THE BULLETIN OF 


Mathematical 
BIOPHYSICS 


THE UNIVERSITY OF CHICAGO PRESS - CHICAGO - ILLINOIS 


VOLUME 12 
1950 


4 = e hes ; Fs 
PUBLISHED MARCH, JUNE, SEPTEMBER AND DECEMBER, 1950 


PRINTED BY THE DENTAN PRINTING COMPANY, COLORADO SPRINGS, COL 


yim 


THE BULLETIN oe. 


Se 
BIOPHYSIC! 


| MARCH 1950 
_ NoTE ON A MECHANISM OF THE DISTRIBUTION OF WEALTH— 
H. D.Landahl. - - - = (= = = = 5 4 ee 
- OUTLINE OF A PROBABILISTIC APPROACH TO ANIMAL SOCIOLOGY: 
Tii—Anutol Rapoport 2 -" = <a ee en ee ies at ere w 
__ ON THE SINGLE EVENT HYPOTHESIS IN RADIOGENETICS— steer 
I, Opatowski and Alice M. Christiansen - - - - -. - . 19 
A PROBLEM IN RADIOBIOLOGY: DIFFUSION AND RECOMBINATION 
Or IONS —A Gs Landau. 46 = ese a oat eee ee eat a OT 
NOTE ON THE GLYCOLYTIC COEFFICIENT—George Karreman - 35 
A NOTE ON LANDAHL’ § THEORY OF PSYCHOPHYSICAL DISCRIMI- 
-NATION—N. Rashevsky - - - - = = -:.2. 825.5 89 
4 ON THE REMOVAL OF AIR-BORNE. DROPLETS BY THE HUMAN | ae ; 
RESPIRATORY TRACT: I. THE LUNG—H. D. Landahl - =. -. 48 
~ Tur Strapy STATE KINETICS OF SOME BIOLOGICAL SYSTEMS; IIT. 
- THERMODYNAMIC ASPECTS—John Z. Hearon - - ‘= =» “57 
THE STEADY STATE KINETICS OF SOME BIOLOGICAL SYSTEMS: hig 
- THERMODYNAMIC ASPECTS—John Z. Hearon: .- - -: = 85 
# Book REVIEW—GEORGE Boo.e, The Mathematical Analysis Ke ti, 
eee ee Hupopost, Eh keg a el et OE de bee 


“THE ee OF CHICAGO PRESS - CHICAGO - ILLINOIS 
VOLUME 12 aa ) t. Ni og ° NUMBER 14 


Tone BU Lb Eb LN O F 


MATHEMATICAL BIOPHYSICS 
EFDITED BY Ne. RAS HE VS K 7 


The Bulletin is devoted to publications of research in Mathe- 
matical Biology, as described on the inside back cover. 


THE BULLETIN is published by the University of Chicago at the University 
of Chicago Press, 5750 Ellis Avenue, Chicago 37, Illinois, quarterly, in March, 
June, September, December, The subscription price is $7.50 per year, the price 
of single copies is $2.25. Orders for service of less than a full year will be charged 
at the single-copy rate. [Patrons are requested-to make all remittances payable to 
The University of Chicago Press in postal or express money orders or bank drafts. 


' THE FOLLOWING is an authorized agent: 


For the British Empire, except North America, India, and Australasia; The 
Cambridge University Press, Bentley House, 200 Euston Road, London, N.W. 1. 
Prices of yearly subscriptions and of single copies may be had on application. 


CLAIMS FOR MISSING NUMBERS should be made within the month following the 
regular month of publication. The publishers expect to supply missing numbers © 
free only when losses have been sustained in transit, and when the reserve stock — 
will permit. 


RUSINESs CORRESPONDENCE should be addressed to The University of Chicago — 
Press, Chicago 37, Ill. 


COMMUNICATIONS FOR THE EDITOR and manuscripts should be addressed to N. é 
Rashevsky, Editorial Office of The Bulletin of Mathematical Biophysics, 5741 © 
Drexel Avenue, Chicago 387, Ill. 


IMPORTANT ANNOUNCEMENT 


The increase in the output of papers in the field of mathematical biology — 
makes it difficult to insure prompt publication without an increase in the size of — 
the journal. Therefore, The Bulletin of Mathematical Biophysics atee eee the 
following service: } 

Upon acceptance of a paper, the Editor will, if necessary, ask the author to © 
shorten the paper to an extent dictated by the requirements of a reasonably 
prompt publication, The shortening should in no case reduce the paper to a mere © 
abstract. Such a shortened paper will be published within six months or less. — 

The unabbreviated original manuscript will be kept on file at the editorial 
office, Any person desiring to avail himself of the complete manuscipt may ob- 
tain promptly a microfilm copy of the latter, at the cost of 5¢ per page plus a — 
carrying charge of 25¢ and postage, by applying | to the Editorial Office, 5741 
Drexel Avenue, Chicago 87, Illinois. 

All papers in The Bulletin which have been thus shortened will be marked at 
the end. by the symbol MF followed by a figure indicating the number of double- 
spaced typewritten pages of the unabbreviated manuscript. © 


en nn ee ttt 
PRINTED BY THE DENTAN PRINTING COMPANY . .°. COLORADO SPRINGS, COLORADO 


BULLETIN OF 
MATHEMATICAL BIOPHYSICS 
VOLUME 12, 1950 


NOTE ON A MECHANISM OF THE 
DISTRIBUTION OF WEALTH 


H. D. LANDAHL 
COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO 


ees population of individuals is considered to be characterized by a 
distribution of wealth which is the result of gain and loss from the 
outside as well as a result of exchange. The population is assumed to 
be uniform except for the variation in a parameter determining the 
manner of partition during exchange. The parameter is assumed to be 
normally distributed over the population. The resulting distribution of 
wealth is calculated for some special cases and it is found to be very 
asymmetric. 


The mechanism of the distribution of wealth has been repre- 
sented in terms of certain models (Marschak, 1943; Kalecki, 1945; 
cf. also Davis, 1948). The present representation is suggested from 
considerations of interacting neural nets (Landahl, 1947) which can 
account for some types of motivational interaction. However, rather 
than attempt to work out a detailed kinetics of interaction in these 
terms, we shall introduce simplifying assumptions. 

We consider a large population of individuals who are all equal 
to each other except in their ability to gain or lose from an inter- 
action. Let the individuals interact in pairs only and in such a man- 
ner that a fraction of the wealth of each is used in common during 
the interaction of mean duration 7. At the end of the time 7; that 
part held in common is redivided in a manner which depends upon 
the magnitudes x; and x; of the abilities of each individual. We con- 
sider the contacts to be made at random in the population. In each 
case let equal amounts be contributed so that the amount is propor- 
tional to the wealth of the individual having the least. Let the divi- 
sion, at the end of the time 7, be carried out in such a manner that 
the gain of one at the expense of the other is in proportion to the dif- 
ference between the values of the parameters x; and x;. Let W; be 
the wealth of individual 7 who has a parameter x;. Then when 7 in- 
teracts with 7, such that W; < W;, the amount involved is 26 Wi, 
6 being a constant fraction. The loss to individual 7 or gain to indi- 
vidual j is 0 W: y(x; — xi) where y is a constant which determines 
the effectiveness of a given difference in x. It is likely that a para- 
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meter of this type x will to a first approximation be distributed nor- 
mally over the population. Thus if we define p(x) by 


p(x) = 


V2a 
we may_take p(x) to be the distribution of x in the population. The 
units of « can be considered to have been absorbed into the quan- 
tity y since only y x is ever required. 

It should be noted that this linear approximation makes it mathe- 
matically possible for the loss to individual 7 to exceed his wealth. 
In order that this should not occur we must-require 6 y(x; — #1) <1. 
However, if 6» > 0.2, the probability of this occurring would be 
less than one in several millions. 

We shall only consider the case in which the interaction is a ran- 
dom process. Thus the net loss to individual 7 during one contact is 
on the average 


eis, | p(x)dx =1, (1) 


| 87 Wilaj— as) (@)day;, (2) 
while his net gain, from individuals for which x; > 2;, is 
if 6» W; ("41 — %;) p(x;) dx; . (3) 


Let the rate of production of wealth by each individual be 
P+ AW(«x) + Bx so that the production may be correlated with the 
wealth as well as the parameter x. We shall not consider anything 
beyond a linear approximation so that P, B and A are constants and 
p> 0. Let the rate of loss of wealth be given by L. + aW (x), where 
L represents a minimum requirement, and a is a positive constant. 
In an interval 7 the net wealth of an individual changes from W (2 , t) 
to W(«,¢t+ 7) so that the gain is W( x,t +7) — W(x, t) or ap- 
proximately 7dW (xz, t)/et. To this extent we may then write (re- 
placing x; by x and x; by x’) 


a I (2 — 2”) Wa’, t)p (2) de’ 
by [ @—2) WO, typ (ay da’ (4) 


+ Pr+ArW (x,t) + Bra —Lr—arW (a, t) = W(a,t). 
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If we differentiate both sides of equation (4) with respect to x 
with ¢ constant, we obtain 


i W (x', t) p(x’) da’ + W(x, t)z(2) 


a Wo(25t) [ee(r)—p(2) —a] (5) 
“are ate CEE, 
ox ot 
in which 
7 =7/6y; (6) 
a= (a—A)r/Oy; (7) 
z2(x2) = [ p@yae’, (8) 


and W' (x , t) denotes the derivative with respect to x. In obtain- 
ing equation (5), we have used the identity »(2) equals the integral 
of x’ times p(x’) over the range x to infinity. : 
On differentiating again with respect to x , we obtain from equa- 
tion (5) 
ow ow 
— 227) —— er =0. 4) 
oe Ox ‘ aol (9) 
If we set W(x, t) =U(x) V(t), then on substituting into equa- 
tion (9) we obtain, on rearranging, 


{laser pe) —ax2 (2) | 0 — 22 (a) 0/0" —— V/V i 10) 


in which U’ and V’ are derivatives with respect to x and ¢ respec- 
tively. Since the first term is independent of ¢ and the second is in- 
dependent of x, they both must be equal to a constant which is here 
taken to be k. The equation with the first member of expression 
(10) set equal to & can be directly integrated since U is absent. The 
equation from the second member leads to an exponential. Thus the 
general solution for W(«, t) can be written as sum of all such solu- 
tions U(x) V(t), or, if M, and N; are arbitrary constants depending 
only on k, 


[a + p(x) 


W(x, t) 
(11) 


22(C) do 
le EPS tl NS kt /r’ 
=z| arm wo [ Be ep (Cy eh > 
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We shall now consider only the steady state. The solution can 
then be obtained by setting k = 0 in equation (11). This can be more 
directly obtained from equation (10). From equation (11) it follows 
that M, = W(0, ©) andN, =W’'(0, o). Define Ga(x) as the inte- 
gral which multiplies N, in equation (11). Note that G(0) = 0, 
G'(0) = 1. Thus for the steady state solution we have 


W(a,0)=W(0, 0) + W'(0, 0)Ga(%). (12) 


The solution of the differential equation must also satisfy equa- 
tions (5) and (4). For the steady state the left-hand side of equa- 
tion (5) vanishes for any z, so in particular for x =—0. Thus setting 
x = 0 and substituting W(2, o) from equation (12) into (5) we 
find 


W'(0, c) ii Ga(t)p (t)dt + W(0, «) 


(13) 
— (.899 +a)W'(0, 0) +Br=0. 


Similarly we obtain from equation (4) a second boundary condition: 


W'0,0) f & Ga(é)p(é)dé — (P—L)1’ 


+aW(0,0)=0. oe 
Since Ga(&) depends only upon a for a given &, those two equations 
may be solved to determine W(0, o) and W’(0, o) in terms of the 
parameters. 
If we multiply through by »(2)dz in equation (4) and integrate 
from —o to ow , then because of the cancelling of the time due to ex- 
change, we obtain 


(P—L) — (a—A)Wr=dwW,/dt, (15) 


where Wy, is the total wealth, so that in this simple case the total 
varies with time according to 


jas 
Walt) We 0) ee ns eee (16) 
es 


It may be noted that the rate of income, J , may be defined as the 
sum of that which is produced and that which is gained by exchange. 
Thus we may write 


_ OW (x,t) 


I(x, t) = 


+L aw (aya) (17) 
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For the steady state the income is a linear function of the wealth. 
100 


CUMULATIVE PER CENT OF WEALTH 


CUMULATIVE PER CENT OF INDIVIDUALS 


FIGURE 1 


To illustrate the results some numerical examples are given in 
graphic form in Figure 1. Since we shall consider only the steady 
state and since the income and wealth are in this case linear with each 
other, we shall refer to income only. The cases for which B = 0 and 
ja 0.010 a OL 05 Lae 1G a = 2 jo — 1, 
L=0 areshown. Since only relative values are used in the figure the 
values of (P — L)7r/6y were not calculated. To illustrate the type of 
distribution which actually occurs, the 1935-36 distribution of in- 
comes (families and single) in the United States is given in the figure 
also. 

The calculations were made by taking W(0, o) as the unit of 
wealth, then determining W(x, o) from equation (12) together with 
the condition (13). From these values, I(”, «) is obtained by add- 
ing L [equation (17)]. The fraction of individuals which have values 
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of « less than x is 1 — z(x). Their income is given by integrating 
I(x, 0) p(x’) dx’ from —w to x. Dividing this by the total income, 
one obtains the fraction of the total income obtained by these individ- 
uals who are a fraction 1— z(a) of the whole population. 

From the results shown in the figure it can be seen that as a be- 
comes large (e.g., a becomes large or y becomes small), the distribu- 
tion becomes more and more uniform. If a becomes very small, almost 
all of the income over and above that part represented by L/a is in 
the hands of a very few, and almost all of the wealth would be in the 
hands of the same few individuals. 

The above results are illustrative only, since the relationships 
have been highly oversimplified. Only paired interactions have been 
considered and the population has been considered to be stationary in- 
stead of changing. No other random processes have been considered. 
In spite of this, the above mechanism is able to represent, to an ap- 
preciable extent, actual distributions of incomes. 
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OUTLINE OF A PROBABILISTIC APPROACH TO 
ANIMAL SOCIOLOGY: III 


ANATOL RAPOPORT 
COMMITTER ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO 


The probabilities of the emergence of the two kinds of social struc- 
ture in a 3-bird fiock (chain and cycle) are deduced under the assump- 
tion of certain biases acting on the social dynamics of the flock. In par- 
ticular a bias against the reversal of peck order and a bias against en- 
counters of individuais of disparate social rank are considered. Like- 
wise a distribution of an “inherent” fighting ability is considered which 
influences the outcomes of encounters... A functional relation is derived 
between the importance of this ability and the initial probability of a 
chain structure. 


In previous papers (Rapoport, 1949a,b), hereafter referred to 
as I and II, we have studied the probability distributions for various 
types of “social structures” based on a certain asymmetric, but not 
necessarily transitive, relation between each pair of individuals (peck 
right). It was pointed out that if peck right depended on a definite 
quantity of some single variable associated with each individual, then 
the relation would necessarily be transitive and the only types of so- 
cial structure would be simple chains of the form 


A BS Ls (1) 


Structures involving cycles (A > B > C > A) have, however, 
been observed, and this led us to formulate a number of probabilistic 
approaches to the theory of the origin of structures, where the out- 
comes of encounters and thus the emerging of structures is, in part 
at least, determined by chance events. 

In I and II we have examined the distributions of structures 
emerging from the most “unbiased” kinds of encounters where the 
probability of victory is always 1/2 for each individual in any encoun- 
ter. Certainly this is a highly idealized situation. Actually, as has 
been shown (Collias, 1943) many factors tend to favor victory. Some 
of these are intrinsic (characteristic of the individual over long pe- 
riods of time, for example, size of comb) ; some are accidental (pres- 
ent over short periods, for example, moulting) ; some are “histori- 
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cal” (connected with the history of the individual, for example, so- 
cial rank). 

In our approach we are interested not so much in the outcomes 
of encounters between single individuals as in the kinds of social 
structures that are likely to arise in small flocks. An experimental 
situation for this approach would be the study of a very large num- 
ber of small flocks with the view of determining the frequencies with 
which various types of structure occur initially and in the long run. 
In the present paper, therefore, we shall inquire how various types 
of bias influence the probability distribution of those structures. 


Bias against Reversal of Peck Right. 


It was assumed in II that from time to time encounters take 
place between individuals in a flock and that as a result of such en- 
counters the peck right may be preserved or reversed with equal prob- 
ability. Under these assumptions, the ultimate, steady state distri- 
bution of structures was computed. We shall now introduce a bias 
against the reversal of peck right, that is, victory is more probable 
for the dominant individual in any pair. Let us see how this bias 
affects the ultimate distribution. 

The fundamental quantities used in the computation of ultimate 
distributions are the “social mutation” probabilities, a;;, that is, the 
probabilities that as a result of some single encounter, structure S; 
shall mutate to structure S;. Evidently 

1 (4,3) 
yj = D6, PD, (2) 
k=1 
where e;“’?) are the probabilities of those encounters which may affect 
the mutation S;;, and 7, are the corresponding probabilities of peck 
right reversal in those encounters. 

In a completely unbiased situation, all the e;,7) are equal to 
[3N (N — 1)]-, since all encounters are equally probable; further- 
more all the 7, are equal to}. However /(i, 7) varies with the pair 
(7,7). Thus, for the case of three individuals 1(1,2) = 1, since only 
one encounter can affect the mutation S,.; but 1(2,1) = 3 since every 
encounter in a cyclic flock can affect S.,. Therefore in that case 
Aig = 1/6; de, = 1/2. 

In introducing an anti-reversal bias, we are now modifying the 
weighting factors 7. We shall consider first the case of three individ- 
uals. 

Let 6, < 4 and 6, < 4 be the corresponding biases against the re- 
versal of peck right between two individuals whose social rank differs 
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by 0 and 2 respectively. The case where 6) = 6. = 4 reduces to the case 
studied in I, where the structure is completely determined after the 
initial encounters. Therefore we shall suppose that 6.5. < 1/4515 
not both quantities are equal to 1/2. The corresponding bias for indi- 
viduals differing in social rank by 1 is irrelevant, because encounters 
between these individuals do not affect mutations. Then 


Qh: 15/6 + 62/3; (3) 
Qoz = 1/2 — 6y. 
We have seen in JI that the ultimate probability for structure S, 
(simple chain) is given by 


Qe1 
Sate eee 
Therefore 
c= 1/2—3, eg t/2= bs 
ant G ee, 8 


(5) 


= ae ay oe 
This quantity will be greater than 3/4 (the unbiased ultimate prob- 
ability of S,) provided 
: 1—2 6 


Se 2 

2— 3 6) — 62 / 

2—46, > 2—3 bd) — 62; (6) 
Soe Oe 


The result is quite obvious intuitively, since the condition 6, > do 
means that the social ranks of the individuals with 2 and 0 dominant 
peck rights tend to be preserved more than the social ranks of indi- 
viduals with 1 dominant peck right. But individuals with 2 and 0 
peck right are found only in the simple chain structure. Hence the 
simple chain structure tends to be preserved. Incidentally, the con- 
dition that the ultimate structure shall be certainly a chain is 
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IG 3G Be hans 


al bd 


(7) 
3—66,=—4—6 6) —2 63; 


6,=1/2. 


In other words, given the bias here considered a necessary and suf- 
ficient condition that the ultimate structure of the three individual 
society shall be certainly a chain is that no reversal of peck right 
shall ever take place between two individuals whose social rank dif- 
fers by 2. 

It is easy to see that the analogous condition for the ultimate cer- 
tainty of the cycle is 6) = 1/2. The case 6, = 6. = 1/2 is indetermi- 
’ nate. But, as has been pointed out, this reduces to the case treated 
ins. 


Bias against Encounters between Individuals of Disparate Rank. 


The only logical bias in the case of three individuals is against 
the encounter between the individuals with 2 and 0 dominant peck 
right respectively. Then, if s denotes this’ bias, (« < 1/3) we have 


Qi = (1/3 — e)1/2—=1/6— ¢/2; 
M1 = 1/2; (8) 
172 14 3 


= = 8. <n —, 
1/6 — 6/2 + 1/2 Seo AS oe 


Again we see that ultimate probability for S, approaches unity as « 
approaches 1/8, i.e., as encounter between individuals of disparate 
rank becomes impossible. 

If both kinds of biases are operating, the following formula for 
S, is easily derived by completely analogous methods 


1 aa Oo 
4—66,—3e—26.+6862 
This simple case suggests that the biases operating in favor of © 
the ultimate establishment of chain hierarchies may have to do with 
the lowering of the probability of encounter or the probability of 


peck right reversal between individuals of disparate rank. We have 
seen in I that every structure except the simple chain must contain 


Si 


Si=3 
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individuals of equal social rank, that is, possessing the same number 
of dominant peck rights. Now each encounter between such indi- 
viduals resulting in a reversal of peck right increases the rank of the 
Winner and decreases the rank of the loser. If encounters and re- 
versals are confined to individuals of equal social rank, then conceiv- 
ably the society will gravitate toward a structure in which there are 
no individuals of equal social rank, and that is a chain structure. It 
would be artificial, of course, to impose such a condition on an actual 
flock, since it is highly improbable that a bird is able to compare its 
social rank with that of every other bird. It may be aware of its 
particular peck right relations with every other bird in the flock, but 
can be hardly supposed to be able to count the number of dominant 
peck rights enjoyed by itself and by its potential rival. However it 
is conceivable that the attitude of one bird toward another may 
roughly depend on the social rank of the other quite apart from the 
particular peck right relation between the two. Thus if a bird sees 
another bird pecked often, that is, by many others, it may fight it, 
even though the particular peck right relation between them may be 
in favor of the other bird. On the other hand the leading position 
may tend to be preserved, because no bird ever sees the leader pecked 
at all. All these tendencies favor the ultimate establishment of a 
chain hierarchy as shown mathematically in the simple case of three 


individuals. 


Distribution of an Inherent Quantity. 


Suppose now that each individual is endowed with a certain 
quantity of some property « which is important in fighting. It is by 
no means certain that the individual with the greater x will always 
be victorious in an encounter but his chances are in general the bet- 
ter the greater the excess of his x over that of his opponent. We may 
make the probability of victory of individual X, over that of X, a 
function of (, — #2.) where x, and x, are the respective fighting quali- 
ties of X, and X,. This treatment is quite analogous to H. D. Lan- 
dahl’s treatment of the probability of guessing the stronger of two 
stimuli in a discrimination test. (Landahl, 1938, 1939, 1940.) 

Let now F'(xz) be the distribution of the quantity x in the entire 
population which is to be divided into flocks. To fix ideas we take 
our flocks to be triples. The range of x can be any convenient one, 
and we shall take it as 0 = x = 1 so that F(x)dz is the fraction of 


the total population with x lying in the interval (x, « + dx) and 


[ F@ae=t. 


12 ANIMAL SOCIOLOGY 


Then the probability that three individuals X , Y and Z char- 
acterized by x, y and z respectively have been selected to form some 
flock will be F(a) F (y) F (z) dadydz. Suppose V(x — y) is the prob- 
ability that X wins an encounter with Y. We seek the expression 
that a triple becomes a chain as a result of the three initial encoun- 
ters. A chain can be formed as a result of the following mutually 
exclusive events 


(SY) UA SY) 7 or SY SZ) eZ) a 
9 
or. GX) (Zee 
Here X > Y means X wins over Y. Note that the necessary and suf- 
ficient condition for a chain in a triple is that one individual is pecked 
by the other two, where the peck order of the other two is immaterial. 
Hence the probability of a chain in the initial distribution is given by 


iP (’ [ Via—w Vie— a) F (@) FW) F (2) dadydez 
+ i f fe V(a—2z)V(y—2z)F (2) F (y) F (2) dudydz (10) 


sl [ J ve-ove—oF @r WF @dvavde. 


To solve an actual problem, we must, of course, make assump- 
tions on the forms of F and V. We shall make the simplest possible 
assumption on F’, namely that any value of x is equally likely in a 
sample. F(x) is then a “rectangular” distribution: F(z) = 1, 
0 =x £1. As for V(&), where é is (wx — y), (z — y), or whatever 
the case may be, we would like it to be of sigmoid form such as the 
usual case with error functions. For large positive values of &, V(é) 
should be unity or nearly unity; for large negative values V(é) 
should approach zero, while for small absolute values of €, V(&) 
should be nearly 1/2 , since those are the cases where victory is a mat- 
ter of chance. 


& 

One such function is V(é) = Hi edt, which is expressible in 
terms of the incomplete gamma function. However, this function 
does not lend itself to integration in closed form. We wish, however, 
to solve some cases in closed form in order to illustrate the method. 
We shall therefore devise some rough approximations. 
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V=0 


FIGURE 1 


Polynomial Approximations. 


Our sigmoid function has the general characteristics of Figure 1. 
The range of é is —1 S & = 1, since the ranges of x, y, and z are 
all (0, 1). Let us find a polynomial approximating such a curve with- 
in the range of €. We assume that the victory of the individuals 
with the highest x over that of the individual of the lowest x is a 
certainty. Therefore if f(€) is the desired polynomial, we must have 


f(—1) =0; f(0)=1/2; f(1)=1; 
f(-1)=f(1)=—0; f’(0)=0. 


To satisfy these six conditions, a polynomial of at most sixth degree 
is sufficient. However, the symmetry of the curve with respect to 
the point (0, 1/2) helps to determine it more precisely. By the 
method of undetermined coefficients one can easily derive the poly- 
nomial 


(11) 


Deady eine om 0/4 Seer (1 3/4) eS, (12) 


where m is an arbitrary parameter. This polynomial satisfies the con- 
ditions (11). If we further demand that f(&) have no turning points 
between — 1 and 1 and no inflexions in that range except at zero, 
this puts some restrictions on m. It is sufficient to demand that no 
inflexions except at & = 0 occur between — 1 and 1 since the exist- 
ence of turning points for this polynomial in that range implies the 
existence of inflexion points. Taking the second derivative, we have 


7. (2) —— (12m — 15/2) é = (20m — 15) &. (13) 


One of the roots of (13) is zero and this represents the inflexion at 
& = 0. The other roots are given by 


[ 24m — 15 
See yyaae Pe (14) 
40m — 30 
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and we wish these to be greater than unity in absolute value or else 
imaginary. In the first case, the range of m can be seen to be 
3/4 < m © 15/6; in the second it is 5/8 = m S 5/4. Since both 
ranges of m preserve the shape of the curve in Figure 1, we may al- 
low the total range 5/8 = m = 15/16. 

Evidently m is the slope of the curve at = 0. The larger the 
m , the more abrupt is the change in the probability of victory of X 
over Y as the quantity (x — y) changes sign. Hence the quantity 
m is a sort of measure of the importance of the inherent fighting 
ability in determining victory. If m is large, the results of outcomes 
are more determined; if it is small the results are more random. It 
would be interesting to establish the functional relation between m 
and the probability distribution of the two structures of triples 
after the initial encounters. However, the polynomial approxima- 
tion allows us only a comparatively narrow range for m. For this 
reason, we shall compute the probability distribution only for the 
simplest case, m = 3/4 and take up the discussion of the functional 
relation between m and the probability distribution in connection 
with another approximation. 

If m = 3/4, f(€) = 1/2 + mé — (2m — 5/4) E*. We wish to 
compute the quantity 


[ [J @-ote—w aedyde. (15) 
We perform the integration with respect to z first. 


[ fe—nie-—nu=te-w [ te—wae 


= f (4 — y) [18/16 — 6/2 + 3y2/8 + y3/4] . (16) 


Next we perform the integration with respect to x. We note that 
the bracket in (16) is independent of x, while f(« — y) is exactly 


of the same form as f(z — y). Therefore the result of the two in- 
tegrations is 


18/16 — y/2 + 3y2/8 + y°/4. (17) 


Integrating (17) with respect to y, we obtain after elementary but 
laborious computations the value .285 , which is 1/3 of the probabil- 
ity of a chain. Hence the probability of a chain is in this case 


S, = .855, (18) 
which is significantly larger than the random probability 3/4. 
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Absolute Exponential Approximation. 


We shall now let our sigmoid curve be represented by the follow- 
ing function: 


V(é) = ters for -&=0; 

V(é) =1—4e for: -€2 0. ee 
This function is continuous and has a continuous first derivative but 
a discontinuous second derivative at & = 0. It is seen to be the in- 
tegral of the absolute exponential function e™é|, which has been 
used by H. D. Landahl (1938) as an approximation to a normal dis- 
tribution. The function V(é) is a good approximation to the error 
function. It does not have the values 0 and 1 at € = — 1 and 1 re- 
spectively, but this simply means that we are abandoning our hy- 
pothesis that victory of the strongest over the weakest individual in 
the population is a certainty. This is equivalent to choosing a smaller 
sample. The range of m is here quite unrestricted (0 = m < o). The 
extreme cases are instructive. For m = 0, victory is arbitrary, and 
the case reduces to that treated in I. For m = ow, victory of the 
stronger individual is certain, and the probability of the chain must 
be unity. Thus for m = 0, we must have S, = 3/4; for m= wo, 
S, = 1. These values will serve as a check on our computations. 

We are interested in the quantity 


fr 7 [Ve@—wVe—wacayae, (20) 


which gives 1/3 of the probability of a chain. To compute the inte- 
gral (20), we divide the region of integration into four nonover- 
lapping regions in which the functions V and limits of integration 
are as follows: 


Region V(*a—y) V(z—y) Limits of « and z 
1. (a@—y) <0; (e—y) <0 semen 2em(2-y) OSuSy; 0S25y 
2. (x—y) <0; (¢—y) >0 semen 1—se"@) USeSy; yS2zS1 
8 («—y) >0; (z—y) <0 1—se-me Zemew ySuSl; 0OS2z5y 
4, (x—y) >0; (c—y) >0 1— semen 1— seme ySeSl; yS2z51 


The range of y is 0 = y =1 for all cases. Thus we must evalu- 
ate the following sum of integrals 


16 ANIMAL SOCIOLOGY 


1 Oy, gy 
{ { f 1 eme-emy+me d 4 dzdy 
Q 0 0 


ff a f jem ies deme] dadzdy 
sees (21) 


© ym 1 
+ J J i) [1— fermen] geme dadady 
0 Os gy, 


4 J | J | J ‘[1— geome] [1 — deme] dardady . 
) y y 


The integrations are laborious but straightforward and yield in 
the last result 


9 — 4m — 2e-" (4 + 8m) — e™ 
ee (22) 
Am 


Multiplying expression (22) by 3, we obtain the probability of a 
chain as a function of m, namely, 


27 — 12m — 6e-" (4 + 8m) — 3e2" 
4m> 3 


We note that Lim S, = 1, as should be the case. For m — 0, the 


S,=1+ (23) 


expression is indeterminate, since both the numerator and the de- 
nominator vanish. However it can be easily evaluated by L’H6pital’s 
rule after three differentiations. Thus upon differentiating the nu- 
merator and denominator of the fraction in (23) once we obtain 


—12 + 6e(4 + 3m) — 18e™ + Ge 


. 24 
12m? 9) 
which is still indeterminate for m= 0. 
On second differentiation, we get 
— 6e"(4 + 38m) + 18e + 18e- — 12e-2" 
————————_- (25) 
24m 
still indeterminate. 
But finally on third differentiation we obtain 
6e-™ (4 + 8m) — 54e-™ + 24e2™ 
pe esr ge ee (26) 


24 ; 
whose limit as m approaches zero is —1/4. Thus Lim S, = 3/4, as 


m—->0 


desired. 
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The above method illustrates how a functional relation between 
the importance of inherent characteristics and the distribution of 
social structures can be found. The actual experimental problem is, 
of course, far less determinate. In general what will be observed is 
the distribution of structure types, from which it is desired to deter- 
mine the bias parameters. We have seen how several bias parameters 
may be responsible for the deviation of the actual structure distri- 
bution from the random distribution. Obviously a single experi- 
mental situation will not enable us to say which bias factors are ac- 
tually at work. However if independent observations are made in 
societies of, say, 3, 4 and 5 individuals, the problem becomes much 
more definite, since different assumptions on the biases for the case 
of three individuals may lead to widely different conclusions as to 
what will happen in flocks of 4 and 5 individuals. 

The 4-bird flock already shows marked differences from the 3- 
bird flock. In the latter the probability of a chain is a priori high 
(8/4). In the former the a priori probability of a chain is quite low 
(1/8). Thus if a predominance of chains is observed in 4-bird flocks, 
one must conclude the operation of biases. On the basis of what is 
assumed about the biases operating in a 3-bird flock, certain pre- 
dictions may be made about the dynamics of the 4-bird and the 5- 
bird flocks. Eventually a theory may be developed which would de- 
scribe the actual evolution in time of the structure distribution acted 
upon by various sorts of bias (so far we have considered only initial 
and final distributions). It is hoped that these studies may help in 
the development of a mathematical sociology. 

This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of The University of Chicago. 
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It is shown that the assumption of a cumulative action of many in- 
dependent random events as the primary cause of induction of a muta- 
tion fits the experimental data at least equally as well as the single 
event hypothesis. The theory is illustrated on the existing data of X-ray 
induced lethals in the X-chromosome of Drosophila. Possible shortcom- 
ings of the single event theory in connection with these experimental 
data and alternative ways of its modification are indicated. 


The idea that a gene mutation induced by an ionizing radiation 
is due to a single primary event has its basis mainly, if not exclu- 
sively, in the quantitative relationship between the amount of energy 
used, D , and the detected frequency, P, of the mutation. In fact, if 
the experimental values of P are plotted against D, the result can 
be frequently well fitted by the curve 


P(D) =1—exp(— kD), (1) 


k being a constant. In general the range of kD is so small that the 
graph of P(D) appears as a straight line (cf. e.g. Timoféeff-Res- 
sovsky, 1937, pp. 67, 123; Sturtevant and Beadle, 1940, p. 210; Catche- 
side, 1948, pp. 327-330). Considering the energy D as a very large 
number of ionizing particles or ionizations and interpreting equation 
(1) probabilistically (cf. e.g. Uspensky, 1937, pp. 135-187; Cramer, 
1946, p. 203), we see that one ionizing particle or one ionization is 
sufficient to cause the observed mutation; if the energy D is meas- 
ured as the number of ionizations per unit of volume, & in formula 
(1) is the size of the volume within which the ionization must occur 
to produce the mutation. 

U. Fano and M. Demerec (1944) warned against an unqualified 
acceptance of the concept of a single all-or-none event as a primary 
cause of a radiation induced mutation; the latter may appear as ini- 
tiated by such an event only because the observable genetical effect 
is of this kind by its nature: either there is a mutation or there is 
not, although mutations of the continuously varying or of the quan- 
titative character may perhaps offer an exception to this latter view- 
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point. In a recent paper (Opatowski, 1950) the author has pointed 
out that other interpretations of experimental results equally plau- 
sible, but incompatible with the single event theory, can lead to a 
relationship between P and D, which within the range covered by 
the existing experimental data does not differ from formula (1) by 
more than the experimental errors. The purpose of this note is to 
obtain such an interpretation through a different mathematical pro- 
cedure and to apply it to different experimental results. 

Let us assume that the primary cause of the observed mutation 
is a cumulative action of a very large number of mutually independ- 
ent events, e.g. ionizations. The central limit theorem of the theory 
of probability translated into non-mathematical terms tells us that 
if a phenomenon is a consequence of a very large number of mu- 
tually independent random causes, then the quantitative character- 
istic of this phenomenon, measured e.g. by the probability of occur- 
rence of a certain event, follows the normal law of probability (Cas- 
telnuovo, 1925, vol. I, p. 160; vol. II, p. 187; Uspensky, 1937, chap. 
XIV; Cramér 1946, pp.- 213-220). In more precise terms we can 
state this theorem and apply it here in the following manner (cf. 
Opatowski, 1950): Let us consider a gene in the irradiated sample. 
Let the first event (e.g. the first ionization) be induced in that gene 
* by a dose A,D, the second event by a dose 4,D , the third by A;D , 
etc. Since each of these events either consists of or is due to an ioni- 
zation and since the occurrence of an ionization at a particular point 
in the sample is a random event, each of the 4;D’s is a random vari- 
able insofar as the induction of the ith event in the gene is concerned. 
It is justified to assume that all these random variables 4;D ((=1, 
2,3,---- 7) are mutually independent and that they have the same 
probability distribution; there are no reasons for making the oppo- 
site suppositions. Consequently the distribution of 


> AD 

(aed 
tends to be normal as 7 tends to infinity (central limit theorem, see 
e.g. Cramér, 1946, p. 215). If we apply this theorem here by assum- 
ing that the mentioned events are caused by absorptions of some 
amounts of the administered energy, then the probability of induc- 
ing a mutation in a sample through an increase of the amount of 
energy absorbed from D to D + dD is: 


dP = A exp[— h?(D— M)?] dD, (2) 


where A, h and M are constant. Each value of D gives a certain 
frequency of mutation P. The latter is not zero when D = 0, be- 
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cause of spontaneous mutations. If we indicate with S the dose that 
would have to be given to the sample to induce in it the genetical 
effects existing in it in a natural manner and related to the particu- 
lar mutation considered, if the sample were completely free from 
those effects, then the theoretically significant range of D extends 
from — S up to the maximum value of the dose actually used in the 
experiment. Whereas the frequency of spontaneous mutations is a 
genetical measure of the degree of transformation which the popu- 
lation undergoes, the quantity S gives an idea about the magnitude 
of the energies which may be involved in this process. At the pres- 
ent stage of our knowledge about the mechanism of spontaneous mu- 
tations this or any other interpretation of S must be very inconcrete: 
it is expressed in terms of ionizing energy, which—at least in organ- 
isms of not very long life span—is known to be insignificant for the 
process of spontaneous mutation. However S is needed for the pur- 
pose of a mathematical analysis which is here carried out. 

If we assume that an infinite amount of energy D = o would 
induce mutation in the entire sample, i.e, if we put P = 1 for 
D = o, then we obtain from equation (2) as the probability P(D) 
of inducing mutation during the administration of a dose D 


P(D) = fi exp [—h?(a — M)?] da/ tie exp [—h? (a — M)?] dx, 


or 
1—P(D) =[{1—A(hD—hM))/E , (3) 


where K=1+ H(hS + hM) and H(z) is the error integral (WPA, 
1941), and 


H (2) =(2/Va) [exp (—28)de. (4) 


The above formulae give the probability of spontaneous mutation 
for D = 0. Except for a displacement of the origin of the system 
of coordinates and for different units of length of the two axes, the 
function P(D) given by equation (3) and the integral H(D) are 
identical. 

Considering now the original question as to the correctness of 
the hypotheses upon which equation (1) is based, we may ask to 
what extent do the functions P(D) as defined by equations (1) and 
(3) differ from each other in the range in which P(D) is experimen- 
tally known? The question is not easy to answer in this general form, 
because three parameters h, M and S appear in equation (3). How- 
ever we have noticed already that, except for a linear transforma- 
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FIGURE 1. The upper diagram to the left represents the error integral H (x) 
and the curves 
E (gx) =1— exp (— gz) 


for g = 1 and g = 1.6. Among all the E(g)’s the curve H(1.6a) is about the 
best approximation to H(x) (see the Addendum to this paper). 

The lower diagram to the right gives the frequency of lethal mutations in- 
duced by X-rays in the X-chromosome of Drosophila. The dots represent the data 
of I. Schechtman (19380), W. P. Efroimson (1981), C. P. Oliver (1932) and N. 
W. Timoféeff-Ressovsky (1934). The crosses are weighted averages calculated 
by K. G. Zimmer (1934). Crosses overlapping with dots have been omitted. The 
curve is calculated from equation (8) with S = 42, M = 6600 and h = 0.00008. 


tion of coordinates, the function P(D) of equation (3) is identical 
with H(hD). For instance, in the very particular case in which 
M = S=0, we have P(D) = H(hD). A comparison of the graphs 
of the function H (a) and of the function P(x) as expressed by equa- 
tion (1) for k = 1, shows that they do not differ by very much (see 
Fig. 1). The difference between H (hx) and the function P(x) of equa- 
tion (1) can be further reduced by choosing an appropriate value of 
h for each value of &. An analysis in an Addendum to this note 
shows that over the whole positive range of D the approximation 


H(hD) ~1—exp (—kD) , with k/h=1.6 (5) 


is about the best. For sufficiently small values of D the relation 
H(hD) ~ 2hD/\/x holds, consequently in this range k/h ~ 2/\/a = 
1.13. It is seen from the figure that the difference between the two 
sides of equation (5) is not large even over the whole infinite range 
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of D if suitable parameters k and h are used. This difference can 
be made still smaller over partial ranges of D. This fact suggests 
on mathematical grounds that drawing conclusions as to the type 
of the primary process of a radiation induced mutation from ex- 
perimental data on mutation frequency and dose may be misleading, 
since opposite types of processes, that of a single random event and 
of many independent random events may imply the same quantita- 
tive relationship, within experimental errors. 

To test this possibility in a more concrete manner the figure 
shows a well known set of sex-linked lethals induced in Drosophila 
by X-rays, fitted by equation (3) with S = 42, M = 6600 and h = 
0.00003 , where roentgen is taken as the unit of the dose. The fit is 
about as good as those known in the literature based on equation (1) 
[see e.g. Timoféeff-Ressovsky, 1937, p. 67; Sturtevant and Beadle, 
1940, p. 211]. All of these fits including the present one are very 
close to a linear one, within the range of D up to about 6000 roent- 
gen, and all show about the same magnitude of deviation from the 
experimental data at higher doses. The recently published data of 
W. P. Spencer and C. Stern (1948) on the lethal mutations in Droso- 
phila, which are not included in the figure, can be fitted quite well 
with the same values of S and M, and a value of h = 0.000025 (Opa- 
towski, 1950). From a statistical viewpoint the distribution (2) as 
it is used in the present paper is a truncated normal distribution. If 
it were not truncated at D = — S but used up to D = — o , the con- 
stant h would be inversely proportional to the standard deviation. 
This may explain why the experiments considered in this paper give 
a higher h than the fit of W. P. Spencer’s and C. Stern’s results alone, 
since the latter involve a smaller number of experimental data and 
consequently may show a larger standard deviation. 

It should be noted that there is no function (1) which would fit 
all the data of Figure 1 on the lethal mutation in the X-chromosome 
of Drosophila; the lower curve of Figure 1 is not only a fit by equa- 
tion (3) but is also about the best fit based on equation (1). Conse- 
quently if the mutation frequency of 18.3% at 9000 roentgen found 
by I. Schechtman and W. P. Efroimson is correct, the single event 
theory would not hold for the mutation considered. Two facts have 
to be kept in mind in this connection: 

(a) The dose of 9000 roentgen is very high, since 10,000 roent- 
gen give a complete sterilization according to W. P. Efroimson 
(1931). Above the range of about 6200 roentgen there are no other 
- experimental data on the mutation considered, except that just men- 
tioned of Schechtman and Efroimson. The latter fit well to a smooth 
curve together with the remaining data (Fig. 1), including those of 
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C. P. Oliver at 6160 roentgen and N. W. Timoféeff-Ressovsky at 6000 
roentgen (cf. K. G. Zimmer, 1934). The other data of Efroimson 
and Schechtman, for instance, at 4500 roentgen with 11.5% muta- 
tion frequency agree very well with the figure of Timoféeff-Ressovsky 
of 13.77% at 4800 roentgen (cf. Zimmer, 1934). Consequently there 
are no indications that the results of Efroimson and Schechtman 
might be incorrect. 

(b) The complete set of data as they are available today perti- 
nent to the lethal mutation considered is in disagreement with the 
theory that mutation is induced by a single event of one type. How- 
ever these data do not conflict with the assumption that there are 
two different events, each capable of inducing a lethal independently 
of the other. This assumption would give a relation of the type 


P(D) =1—pe"— (1— p)e”, withO<p<1, (6) 


instead of equation (1), and would maintain the essential character- 
istics of the single event theory, except that it would imply that two 
different mutational processes are actually taking place, although no 
differences between the end products of these mutations are detected. 
It is very likely that equation (6) can be well fitted to the data of 
Figure 1; however, such a fit must involve essentially different val- 
ues of k and h, because otherwise equation (6) would be practically 
identical with equation (1). Since & and h are the volumes within 
which the events discussed must take place to be effective, the ques- 
tion arises whether the order of magnitude of k and h, as obtained 
from such a fit, agrees with the presently accepted figures of this 
volume. We leave this question open for the present. Equation (6) 
which has never been considered for Drosophila, has been recently 
suggested for bacteria, on different grounds though than those here 
indicated (Cavalli and Visconti, 1948, 1949). 

The present paper presents some doubts, on a theoretical basis, 
about the validity of the single event theory in radiogenetics and 
indicates the type of mathematical work that still could be done to 
investigate these doubts further. Doubts on the same subject have 
been expressed also on more direct experimental grounds (Caspari 
and Stern, 1948), although they do not seem to be entirely confirmed 
by the more recent investigations of D. E. Uphoff and C. Stern 
(1949). 

This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of The University of Chicago. 


ADDENDUM 
Putting hD=a,k/h=g, E(gx) =1— exp (— gx), R(g, x) 
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ae H(x) — E(g, x), the problem of obtaining the best approxima- 
tion of the type (5) may be stated as the problem of finding a value 
of g such that the maximum value that |R(g, «)| takes as x varies 
be as small as possible. The solution of this problem is dependent 
on the range of x considered. The range 0 to + o is here discussed. 
In order to carry out this discussion the intersection POInts .% —%% 
of the curves H(w) and E'(gx) must be examined first (cf. the fig- 
ure). These points are solutions of the equation: 


gs3—211n [1—A(2,)]. 
The quantity g(«.) is an increasing function of x for x > 0. Its 
graph is very close to a straight line in the range of x from 0 to 2.8. 
It passes through the point x = 0, g = 2/\/a; consequently the 
curves H(z) and E'(gx) do not intersect each other in the positive 
range of x if g <2/Va ~ 1.13. 

We can discuss the problem of the best approximation E'(g, x) 
= H(z) by splitting the range of g into two intervals 0 < g < 2/Va 
and g > 2/\/n. (Negative g’s do not have any meaning in the genet- 
ical problem considered.) Consider first 0 < g < 2/.\/a. In this in- 
terval R(g , x) is positive, and for any constant g the maximum val- 
ue of R(g, x) is attained when z is the positive root x, of the equa- 
tion 

Lm? — JXm +1n(g Va/2) =0. (7) 


Since here 0 < g < 2/\/z, this equation has one and only one posi- 
tive root. The value of R(x, , g) plotted against g appears as a de- 
creasing function of g with its smallest value ~ 0.165 at the end of 


this range, i.e., for g = 2/\/a (the corresponding value of 2%» is also 


BN aL) - 

The approximation just obtained, through g = 2/ Vn, can be 
still improved if g > 2/\/z is considered. In fact, in this range equa- 
tion (7) has two positive roots %» which correspond each to a maxi- 
mum value of |R(x,g)| as x varies and g is fixed. As g is increased 
(beyond 2/\/z) one of these maximum values of |R (x, g)| decreases, 


starting from ~ 0.165 at g = 2/ \/a; the other maximum value, how- 
ever, increases, starting from 0 at the same value of g. Consequent- 


ly under the condition g > 2/ \/n, the best approximation H(x) ~ 
E(g, x) is obtained when these two maxima are equal; this occurs 
about at g ~ 1.6 and gives |R(am, g)| ~ 0.05 (for &m = 1.34 and 
0.26). Since 0.05 is smaller than 0.165 we conclude that g ~ 1.6 gives 
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the best approximation of the error integral H(x) by the function 
E(g, x) =1—exp (—gzx) over the whole positive range of x. 
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A PROBLEM IN RADIOBIOLOGY: 
DIFFUSION AND RECOMBINATION OF IONS 


H. G. LANDAU 
COMMITTEE ON MATHEMATICAL BIGLOGY 
THE UNIVERSITY OF CHICAGO 


_, A perturbation treatment is given for the cylindrically symmetrical 
distribution | of ions or their products resulting from the passage of a 
Single ionizing particle, and taking into account ion diffusion and re- 
combination. The results apply to cases in which diffusion is more im- 
portant than recombination. 


We are concerned with the ion distribution resulting from the 
passage of a single ionizing particle or ray through a homogeneous 
medium. It is assumed that the ionizing particle travels in a straight 
line, and that in the absence of an external field the ion distribution 
is cylindrically symmetrical. The ions or their products* resulting 
from the radiation diffuse and simultaneously recombine at a rate 
proportional to the product of their concentrations. 

The partial differential equation describing this process is non- 
linear so that an explicit solution is hardly to be expected. An ap- 
proximate treatment, that is still used, was first given by G. Jaffé 
(1913) ; for another approximation see A. N. Gerritsen (1948). In 
this paper we treat the problem by a perturbation method, and cal- 
culate the first few perturbation terms, which makes a more accu- 
rate determination of the distribution possible in cases in which dif- 
fusion is more important than recombination. The applications to 
radiobiology and the delimitation of the situations in which these 
results can be applied will not be discussed here, but will be consid- 
ered in later papers by I. Opatowski. 

Assuming that the diffusion coefficient, D, is the same for both 
positive and negative ions or their products, then the partial differ- 
ential equation for ion concentration is 

oP =Dvn—an, CET oO; (1) 


*H atoms and OH radicals in the case of water. 
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where n= concentration of either positive or negative ions, 


t= time, 


lave 
Ves + ——, Laplacian, 
On iP Oa 


r = distance from the path of the ionizing particle, 
a= recombination coefficient. 


For the initial condition, we take, as is usually assumed (Jaffé, 1913; 
Kara-Michailova and Lea, 1940), 


N 5 aed nies 
n(r, 0) ar aa (2) 
Here, N, is the number of ionizations produced by the ionizing par- 
ticle per unit length of path, and b is a measure of the initial radius 
of the ion column. 

In the following we use the method of sources and sinks or the 
method of Green’s function (Carslaw and Jaeger, 1947). This meth- 
od allows us to give an interpretation to the initial condition (2). 
The distribution in (2) is seen to be the same as that resulting from 
from an instantaneous line source of strength N,, acting at time?, 


b2 


in the case of diffusion without recombination. That is, (2) is equiva- 
lent to assuming that the ionizing particle instantaneously produces No 
ions per unit of path length, and that these ions diffuse outward with- 
out recombination for a time, b?/4D , before recombination begins. 
Of course, the actual process is more complicated, involving colli- 
sions off the path of the ionizing particle. 

In order to simplify the expressions which follow it is conveni- 
ent to introduce the following dimensionless variables: 


mb? (3) 
——_— N, 
y N, 
p10, (4) 
D 
Cae (5) 
aN 


(6) 
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Then equations (1) and (2) become 
Ov 
—=V,?»— 2, O<eP< oo; (7) 
Oct 
y(p,0) =e, (8) 


To solve this system by the perturbation method we assume that 
y may be expressed as a power series in f, 


y=) B® ,(p, 7). (9) 
k=0 
Substituting in (7) and (8) and equating equal powers of $, we find 
0 Vo 
= Vp" %, p20; (10) 
OT 
¥(p,0) =e; (11) 
and fork >0, 
0 VK k-1 
= V>)? ¥%,— > Vi Vk-1-j 5 p 2 0; (12) 
OT j=0 
%(p,0) 50. (13) 


Now it is known that the distribution due to an instantaneous 
cylindrical surface source at time, r = 0, which liberates one ion per 
unit area of cylinder and has radius, o , is (Carslaw and Jaeger, 1947) 


Gt 57) =e |—2 | 1 (£5) (14) 
: 27 27 


Ar 


where J, is the Bessel function of zero order and imaginary argu- 
ment. This can be used to write the solutions of systems (10) and 
(11) or (12) and (18) in the form of integrals. The solution of 


3 
a= Vere + yl), o>0; (15) 


y(p,0) =¢(p); (16) 
is 
rp. r)=[ $l) Gl(p,cir)de 
T co (17) 
+ i di [ vlo,A)G(p,asr— Ade. 
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Using this formula the value of » is 
ro(psr) = | €*G(p, 057) do. (18) 
¢0 


To evaluate the integral in (18) we need the formula 


a? 


*c0 i 3 
| ee" Ti(az)edz=—e~ , (19) 
0 ae 


which is readily established by using the power series for J, . 
Then we obtain 


eons (20) 


AUIS Mersey 


This, of course, gives the distribution for o = 0, no recombination. 
It could have been anticipated from the interpretation of the initial 
condition as the distribution due to an instantaneous source, since 
(20) is the distribution due to an instantaneous line source at time, 
+ =-— 1, for diffusion without recombination. It is the first approxi- 
mation given by G. Jaffé (1918). 

From (17) we can now calculate y, as 


¥4(p,7) =Gl— 7 (p, 7) ] = 


ée -_— 
—_ di ——_——_—— G(p,o3;7— A) di, 21 
i I a bey (21) 


where we use G to denote the second integral operator in (17), that is, 
Glw(e.1= | di [ w(e,a)G(p.esr—Aada. (22) 


It is convenient to introduce a new time variable 
C=4,-+1= (4Dt + b?) /b?, (23) 


and correspondingly 


p=424+1; (24) 
then using (19), we obtain after some manipulation 
il E | 2p? d u 
v,=—- exp | ——— ——_———, 25 
pte 2o—wd u(2o—p) a) 


which can be reduced to 
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p?/e 
e —p7/¢ dx 
Lal cae er a ex 
8f ? x 
p?/g(2g-1) 


€ 


> dt 
— Hi(— 2x) =| ere (27) 


is the exponential integral which has been extensively tabulated 
(Jahnke and Emde, 1943; Mathematical Tables Project, 1940). It may 
be readily verified directly that (26) satisfies the equation for »,. 

We have also calculated the next term, »,, which has a rather 
complicated expression. Starting with 


¥,=G[—-2%%], (28) 


(26) 


where 


and using the formula 


foe) box? 
dy re 
e¥— je" I, (ax) x dx 
0 b32? y 
4c2 4c? C2+b2 
e i dz 
= —— (Be 4 
ve 
2¢ a2 by ov 


402 c24b4 
which is readily established by a change of variable and interchange 
of order of integration, we obtain after some manipulation 


Q/e U 
e ue (" dy da 
Y= — ev — | e*—, (30) 
86 p2/e(2g-1) | 4 L 4 2 


(29) 


with 
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ia 
Bese | 


One of the quantities of principal interest is the total number of 
ions per-unit path length, 


N(t)= [ 2arndr=2Ne | rane (31) 


Let the fractional number of ions be 
M(t) =N(t)/No- (32) 


The contribution to M(t) from » and y, could be calculated directly, 
but the use of the following formula is simpler and also makes the 
calculation of the contribution of ». possible. Referring to the sys- 
tem (15) and (16), we obtain by integrating over the region bounded 
by p=0,p=0,7—0,7r=7, 


Jrede= J sdedet+ [oar ['vle.rpdp, (88) 


which is a form of Green’s theorem for this system; and physically 
simply states that the total number of ions is equal to the initial num- 
ber plus the additional number generated by the sources. Then we 
obtain 


My=2 [vp dp=1; (34) 


re) T ee) 2p2 
é PA 
M,(r) =2 v¥,padp——2 di ——_—— = 
1 J pWwp (| f (44+ 1)2 pap 


—tflog(1 +47); 


ee, (oo) _ 29? p2/m x 

lu dx 

mco=a | al e "| f olay 
ei M 0 p2/m(2u-1) v 


Integrating by parts [or using limiting case of (29), a > 0], we find 


(35) 
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1 1447 4u—1 du 
sine gottoen!°S is 


3 ll 
1 i 4 1 du iL. 4 . 
—— log - u + lo (1-= )J*= 2 {[ Ss 
32 J, 3 zg ret Ae oe Loe (1+ 47)] (36) 
== [510 ~) feos 
( e3 39 log u = 


3+16T 


44+16T 


The last integral is a tabulated function (Powell, 1943). 

For large values of +, M.(r) — [4 log (1 + 47)]%, so that 
M = M, + $M, + £?M, approaches the first three terms of the ex- 
pansion of the approximation arrived at by Jaffé (1913), 


a No 4D a 
N/No=| 1+ ] (1+= )|. 
/No 82D og b2 (37) 
D. E. Lea (1947) considers the case in which the two types of 
ion products, H atoms and OH radicals, have different diffusion co- 
efficients, D, and D,, and the columns have correspondingly differ- 
ent initial radii, b, and b.. For this case it is not difficult to calcu- 
late the first two terms of the perturbation calculation of N, by the 
same method as was just used. The result is 
a No D,+ De. 
N/N. = 1 — ————— lo jt 4( |e | 38 

/No aD DS) g b,? + ba? ’ (38) 

which is the same as the first two terms in the expansion of expres- 


sion given by Lea. 
2 


b 
Our previous simplified interpretation of aD as the diffusion 


time before the start of recombination can be applied to the numeri- 
cal values given by D. E. Lea (1947). These give 8 X 10-° sec. for H 
atoms, and 0.5 X 10-° sec. for OF radicals. 


The author is indebted to Dr. I. Opatowski for calling his atten- 
tion to this problem and for helpful discussions. 


This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of The University of Chicago. 
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The glycolytic coefficient of a cell is expressed as a function of the 
external concentrations of glucose, lactic acid and oxygen. It is shown 
that this function has a smallest value when the external concentration 
of glucose is zero and a certain inequality between the external concen- 
trations of lactic acid and of oxygen is satisfied. 


The glycolytic coefficient 6 of a cell is an important constant, 
occurring in the theory of cell division and cancer research as de- 
veloped by N. Rashevsky (1948). This theory raises a number of 
problems which might be attacked by minimizing the glycolytic coef- 
ficient 6. The latter will be the purpose of the present work. From 
the theory of cell respiration of N. Rashevsky (1936) and H. D. Lan- 
dahl (1937, 1939) the glycolytic coefficient 6 can be expressed as a 
function of the external concentrations ¢ 1 , Cox and Cos of glucose, lac- 
tic acid and oxygen respectively. 

Denoting 


Cop =U, Coo—Y, Cor —24; 
GAy ol 7: 
din 7 
pnfh,=s; 
pfds=t; 
af Achy = U; 
ZAsf2p? Aon = V3 
Z2Asf*pAs’n a=W; 


[in which a, f, A:(i=1, 2,3) have the same meaning as in the first 
above-mentioned work], we find as the result of elementary calcula- 
tions: 

p 
B Se SS 
3(py + 72) 


(1) 


| 2rz ~y| 
sx —ty —uz—p + VV (sx + ty + uz + p)? — 2vxry — 2waz 
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Denoting further for abbreviation: 
settyt+tuet+pHA(z,y,2)3 (2) 
V (sa + ty + uz + p)? — 2vxy — 2waz= B(x ,y, 2); (3) 


$i. — ty —uUz — @ 


(4) 
+ V (sa + ty + uz + p)? — 2vay — 2waxz=C (x,y, 2); 
we get 
2, 
p=? (= -y): (5) 
3(py'+7z) \ C 
Simplification of 
0 e 0p 
50; “F =o; apigaee (6) 
Ox oy i 
leads to respectively: 
A — vy — wz 
B 
tA — vx 
9 2rz (- pS S| 
Zz 
ee (—=-v)- SS eee ee (8) 
py t+ rz G C2 
UA — wx 
9 ar —2r2(—w+=—"* 
r 
——?_( 2% _,)+____-___8 "=o. (9) 
py + rz’ C C2 


Unfortunately, it is impossible to satisfy these three equations 
(7), (8) and (9) simultaneously. However, it is easy to see that 
z = 0 satisfies the first and second equation, while the rather wide 
condition sz — ty > 0 or, what amounts to the same thing, 


assures that the left-hand member of equation (9) is positive. There- 
fore, in this case 


0p 0B 0p 
——03. —-—0,— -——. 0 


Ox Oy 02 


so that 6 has a smallest value when z is as small as possible. The 
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latter then has the value z = 0, because the physical meaning of z 
allows it only to have values which are at least 0. 


As 


Concluding, we see that when z = 0 and «3 > Coo, P has a 


J 


smallest value, being — 4. 

The meaning of the negative sign of the glycolytic coefficient has 
been discussed by N. Rashevsky (1948, p. 74). 

This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of The University of Chicago. 
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In his theory of psychophysical discrimination H. D. Landahl as- 
sumes that because of the symmetry of the neural circuit considered the 
fluctuation of the difference of excitations at the two parallel connec- 
tions is the same as if the fluctuations occurred only at either one of the 
connections. It is shown that if the excitation at both connections fluc- 
tuates independently according to a simple exponential probability fune- 
tion, the difference of excitation fluctuates according to a different prob- 
ability function. 


In his theory of psychophysical discrimination H. D. Landahl 
(1937) [cf. also Rashevsky, 1948] assumes random fluctuation of 
excitation in the two parallel cross-inhibited pathways. Because of 
the symmetry of the whole structure a positive fluctuation on one 
side is equivalent to a negative fluctuation on the other side. There- 
fore Landahl considered first fluctuations at only one of the two 
parallel connections. 

In a subsequent paper (Landahl, 1940) it was shown that if the 
fluctuations are distributed according to the Gaussian law, then the 
results are the same, regardless of whether one set of fluctuations 
occurs at one connection or two independent fluctuations occur in 
both parallel chains. 

Inasmuch, however, as a distribution different from the normal 
has been introduced by H. D. Landahl, and found more convenient 
for deriving closed expressions, it is interesting to investigate the case 
of two independent fluctuations for such a distribution. 

The purpose of this note is to show that in this case the restric- 
tion to one fluctuation is not always justified, although it does not 
necessarily invalidate Landahl’s results. In fact what is essential in 
Landahl’s theory is to assume that the difference, 


= €; — &2; (1) 
fluctuates according to some distribution function p(~), such that 
39 
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(e p(x)dx , (2) 


and where p(x)dx denotes the probability of a fluctuation between 
x and z+dz. This does not mean, however, as we shall see, that 
either «, or €, individually fluctuates according to the same function 
p(x). 

Let 4 > 0 for the sake of definiteness and consider the probabil- 
ity P, of the (wrong) reaction R. when both e, and e, fluctuate in- 
dependently according to some function p(x). 

Suppose that at a given moment «, has a fluctuation ¢, while é. 
has a fluctuation ¢’, so that instead of «; — e. we now have (e, + ¢) — 
(ec. + ¢’). Considering for simplicity the case of two categories, that 
is, h = 0 (Landahl, 1937; Rashevsky, 1948), we find that for R, to 
occur we must have (e, + ¢) — (& + ¢’) < 0, or 


ope A. (3) 


Hence for a given ¢, any value of ¢’ which satisfies (3) will result 
in the reaction R,. 
But, for a fixed ¢ , the probability of inequality (3) is given by 


” (ade =F(A + $). (4) 


A+ 


Since the probability of ¢ being between ¢ and ¢ + d¢ is p(¢)d¢, 
the probability of R, for a given ¢ is equal to 


F(A + $)p($)d¢. (5) 


Since ¢ may have any value from — o to+o, 


Po= [ F(A+4)0(y)dg= [ v(g)dg [" p@ar. ©) 
-0 -0 A+ 


Let us evaluate P, under the assumption made by Landahl that 


k 
p(x) a ers (7) 
Then we have 
k 
forx>0 p(x). +p, (a) ti: en (8) 
k 
fora <0 p(x) =p_(x) Feige er. (9) 


For d+ ¢>0, or for ¢ > — A, we find from (4): 
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ice} 


k 
F(z +4) =F.(A+ ¢) =| eda =f eras) (10) 


A+o 


For 4+ ¢<0or¢ < —A we find 


rasy=ra+y a4 | ears 5 | ede 
A+o 2 0 
(11) 


=1— fF et), 
Remembering (8), (9), (10) and (11) we see that (6) now be- 
comes 
eal 20 
Po= | E(A+$)p-()dg+ [P(A + 9)9-(a)de 

-2 -A 

é. (12) 

all F(A + $)p.(¢)d¢. 


Introducing (8), (9), (10) and (11) into (12) we find after 
rearrangements 


1 kA 
Po=( aia ri ) es. (13) 


The quantity P, decreases monotonically with A since 


ra een TA <0) (14) 
cin ee ‘ ; 

Expression (13) differs, however, from that derived by Landahl, 
which is 


Py, =e, (15) 


Thus we see that if A as a whole fluctuates according to the law 
given by (7), the probability P,, is less than if ¢, and e fluctuate in- 
dependently according to (7). 

After reading this manuscript Professor H. D. Landahl has 
shown by a simple graphical construction that equations (13) and 
(15) can represent a given empirical curve equally well with a maxi- 
mum difference of approximately 2%. 

The author is indebted to Dr. Landahl for reading the manuscript 
and for the above contribution. 
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Calculations as to the fate of air-borne particles in the respiratory 
tract are given in a paper by W. Findeisen (1935). It is the purpose of 
the present paper to extend and modify these calculations. The results 
obtained for total retention in the lungs are found to be in general agree- 
ment with published data (cf. Landahl and Herrmann, 1948). The cal- 
culated results also give information on the amounts deposited in vari- 
ous regions as well as the amounts retained in different fractions of 
expired air. 


Since the geometry of the respiratory passages is very complex, 
it is at best only possible to represent the important features and then 
assign approximate values to the average dimensions from observa- 
tions. Since these values are very approximate it would seem neces- 
sary only to carry out calculations of the effect of the factors involved 
in an approximate manner. Furthermore, since several factors are 
acting to different degrees in various parts of the respiratory tract, 
the net result does not generally depend too much upon the effect of 
one factor of removal in any one region. For this reason we shall 
derive approximate relations for the three principal factors involved 
in the removal of particles by the respiratory tract. 

Schematic Representation of the Respiratory Tract. In order to 
make an estimate for the probability of removal of any particle from 
the air stream, it is necessary to set up a model of the respiratory 
tract. In so doing, we follow the model of W. Findeisen (1935) fairly 
closely. The respiratory tract is represented by a tree in which re- 
gions of the same order are taken to be equivalent. In most cases the 
regions closely approximate cylindrical tubes. The “mouth” and 
“pharynx” are added, and two orders of alveolar ducts are used in- 
stead of one in order to reduce the very large branching factor. A 
few values are slightly changed to conform to measurements by vari- 
ous other investigators. But the principal changes are made to sat- 
isfy the requirement that the volumes calculated for the alveolar 
ducts and sacs be sufficiently large (cf. Willson, 1922; Loosli, 1949). 

«The work described in this paper was done in part under contract between 
the Medical Division, Chemical Corps, U. S. Army, and the University of Chi- 


cago Toxicity Laboratory. Under the terms of the contract the Chemical Corps 
neither restricts nor is responsible for the opinions or conclusions of the author. 
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The numerical values of the geometrical measurements for the suc- 
cessive orders, a, are given in Table I. From these other values are 
calculated for a flow rate of 300 cc/sec., inspiration time of 1.5 sec., 
and 450 cc. tidal air. The last columns will be discussed subsequently. 

Removal by Impaction Within Respiratory Passages. Because of 
its inertia, a particle suspended in an air stream will not be able to 
follow the stream as it bends around in its passage through the res- 
piratory tracts; as a result, it may impact against the moist surface 
of a passage and be removed from the air stream. In this section we 
shall not consider impaction on cylindrical objects such as the nasal 
hairs. In the case of a tube, bent at 90 degrees, it is found (Landahl 
and Herrmann, 1949; cf. also May, 1945) that the fraction of the 
particles removed from the air stream is about one-half when the 
following relationship holds: x = pd? v/18nR = 1, where p is the 
density and d is the diameter (cm) of the particles, v is the average 
velocity of the air stream, 7 is the viscosity of air (1.9 X 10-* g cm 
sec.+), and where R is the radius of the tube. 

If we have a tube which bifurcates or branches, then the impac- 
tion in such a tube of order a, will be about fifty per cent when 
pd?Va sin 6/18yR, = 1. The quantity @ is the angle between the 
axes of the tubes of order a — 1 and a. The quantity v'c, sin 6 is 
the component of the velocity at right angles to the second tube. Ac- 
tually, v’c, is not exactly the same as va, the average velocity in 
a tube of order a — 1, but the average velocity of that part going 
into a tube of order a. Unless the number of branches at one point 
is very large, the discrepancy is not serious, and v’., is taken equal 
to Va. throughout. The average distance that the particle must pass 
perpendicular to the air stream is now R., the radius of the tube a, 
instead of FR. 

From the values given in Table I, it can be seen that the Rey- 
nolds’ number (cf. Lamb, 1932) is in the general neighborhood of 
that which is required for turbulent flow. Thus we would expect the 
actual distribution of velocities to lie between a rectangular and a 
parabolic distribution. If all particles are moving at about the same 
speed, then the fraction which is removed by impacting against the 
outer wall at the bend of the tube would be roughly proportional to 
the ratio of the distance of slips through the air stream to the tube 
diameter. Therefore, in view of the result referred to in the 
preceding paragraph, the proportion would be roughly 3pd?va1 sin 
6/18Ran when this quantity is less than unity. On the other hand, if 
the velocity distribution is parabolic, the fraction which impacts would 
be less than the above linear approximation for very small x or for 
x much greater than unity, but the fraction would be greater than 
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that given by the linear approximation for values of x about equal to 
unity. For example, a particle on the inner wall of the bend has a 
low velocity but must travel twice as far as a particle near the cen- 
ter of the tube, which travels with a higher velocity. These considera- 
tions suggest that if we choose the expression for the fraction re- 
moved by impaction to be of the form x2/(1+x), where x = pd’Va1 
sin 6/187Ra, the result should be adequate, except perhaps for very 
small values of x. When x=1, the fraction is one-half. On the aver- 
age, 6 is about 30°. If Ia‘ denotes the fraction of material in section 
a — 1 impacting in section a, due to inertia, we have 


Py = 150 Pp d? Vaa/ (Ra + 150 p d? Vain) (1) 


as the approximate equation determining the fraction impacted on 
inspiration. On expiration, the amount, I.°, would be the same ex- 
cept for the replacement of va: by va1. In the range of values ac- 
tually used in the calculations, equation (1) gives substantially the 
same values as given by the rather complex equation derived by W. 
Findeisen. 

Removal. by Sedimentation. Due to the action of gravity, a par- 
ticle settles downward through the air. In the absence of turbulence, 
particles of a given size very quickly attain a velocity such that the 
force of gravity is opposed by the force of friction. For particles that 
are not too small, the latter force can be calculated from Stoke’s law 
(cf. Lamb, 1932), but if the particles become small enough to be 
comparable in diameter with the inter-molecular spaces of the air, 
the velocity is greater by a factor depending on the particle radius 
and the properties of the medium. If we write M; for this factor, 
then (Findeisen, 1935) 


Ma=1 +2 X 10° (0.864 + 0.29 e-8-25«10!4) /q , 
or( Millikan, 1947) 


Ma=1+1.8X 10°/d, (2) 
and the velocity of settling is 
: Va = 2.9.X 10° p d? Ma cm/sec. (3) 


vertically downward. Let the time that the particle remains in a re- 
gion a be za. In this time the particle would fall a distance vi ta. If 
a particular tube is inclined at an angle wy with the horizontal, the 
vertical distance that a particle must fall before settling on the inner 
surface of the tube is increased by one over the cosine of the angle. 
The ratio of the distance of a fall, vg 7a , to the maximum distance of 
the fall, 2R/cos y, is then 2.9 X 10° p d? rz Ma cos w/2Rza, so that the 
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fraction which settles out of the air stream will be roughly propor- 
tional to this quantity. 

The average value of cos y for a randomly placed tube in three 
dimensions is 2/4. We shall use this value with a few exceptions to 
be mentioned subsequently. In the absence of any turbulence each 
particle moves through the air a distance of one radius in the direc- 
tion perpendicular to the tube when 2.3 X 10° p d? 71 Mz = Ra. The 
area of intersection of two equal circles, each passing through the 
center of the other, is 39% of the area of either. If each particle 
moved down one radius, then all the particles not in the intersection 
would settle out, that is 61% is removed. Thus, if S is the fraction 
removed in section a due to sedimentation, S = .61 when 0.79 vg ta = 
R.. If one makes the rough approximation that the relationship is 
linear, then a more careful analysis shows that Sz is slightly under- 
estimated for small values of vg ra/Ra, but, of course, overestimated 
for large values of vg ta/Ra. The effect of any turbulence would be 
to increase Sz, slightly for small values of S. but to decrease S, ma- 
terially for larger values. Furthermore, the settling in quiet air, 
which occurs in the region of orders a — 1, leaves the particle dis- 
tribution in the region of order a altered, on the average, in a man- 
ner that slightly favors sedimentation for small S.. If we set 


Sa = 1 — E18 x 105 p @ Ma Ta/Ra : (4) 


then instead of 61% for 0.79 va ta/Ra , we obtain 55%; while for small 
values of Sa we overestimate the settling in a single smooth section 
by a factor of about 1.25. In view of the above discussion, this ex- 
pression for Sa should be a satisfactory approximation. If cos yw is 
not taken equal to 0.79, then 1.8 should be replaced by 2.3 cos y. To 
calculate the sedimentation during inspiration, S.‘, which in the cal- 
culations always equals the value for expiration, the quantity 7a is 
set equal to va‘ (Table 1). For the sedimentation occurring when the 
breath is held, Sa”, 7a is set equal to 7”, the time that the breath is held. 

Removal by Brownian Movements. Small particles suspended in 
air are in continuous motion and after a short time will have moved 
appreciable distances. If D’ is the diffusion coefficient for the par- 
ticle in air, then the mean displacement after a time 7 will be A = 
\/2 D’; (Einstein, 1905). But, since D’ = kT°/3 2d, k being Boltz- 
mann’s constant per molecule, and 7° the absolute temperature (Ein- 
stein, 1905), we have 


A=V 2kT° M./327 d=6.85 X 10° VMaz/d | (5) 


in which the effective viscosity in the neighborhood of the particle is 
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taken to be reduced by the same factor as in the case of settling. We 
may estimate the fraction of the particles removed by diffusion in a 
rough manner as follows: Suppose A = 0.3R. Consider a circle of 
radius R with another concentric circle of radius 0.7R . Assume that 
the particles are moving only in three mutually perpendicular direc- 
tions, one-sixth moving to the left, one-sixth to the right, etc. Con- 
sider the particles on the inner circle and on a horizontal axis. Of 
those moving to the right, one-half will travel 0.8R (since A= 0.3R), 
or 1/6 X 1/2 of the total. Of those moving to the left (1/2)1/* = 0.02 
will be removed, or 0.003 of the total. Of those moving up or down 
(1/2) -7/3 = 0.20 will be removed, or a total of (1/6 + 1/6) X 0.2 = 
0.067. Thus the total removed will be 0.083 + 0.003 + 0.067 = 0.158. 
Since 0.7 is approximately equal to 2', about one-half of all particles 
are initially interior to the circle of radius 0.7R . For those inside the 
inner circle, the probability of removal will be less, whereas for those 
outside, the probability will be greater. Since the mean probability 
is not likely to be much different from the median calculated above, 
we ignore the discrepancy. For not too large values of A/R we can 
expect that the proportion of particles removed from the air stream 
due to Brownian motion‘is roughly proportional to A/R. However, 
for large values of A/R, the probability can approach unity only 
asymptotically. Since the free paths are distributed exponentially, an 
expression of the form 1—e~”, where y is proportional to A/R, is 
likely to be satisfactory. Hence, we may choose y such that for 
P= 0.153, A = 0.3R, and thus obtain for the fraction removed by 
diffusion, Da , in section a: 


Da = 1 — e+ * 10-*vMa Ta/d/Ra_ (6) 


The values of D. for inspiration or expiration, D,‘ , the fraction re- 
tained by diffusion during the pause, will be obtained from equation 
(6) with ra = 1a‘, while D.” is obtained by setting ra =7". The values 
given by equation (6) are substantially the same as those given by 
calculations from the equations of W. Findeisen. 

Wall Effect and Size Effect. Because of its size, a particle cannot 
have its center too close to the wall. This effect can be shown to play 
a very small role. Where the effect might be appreciable, other fac- 
tors will have’generally removed the particle. Similarly, for exam- 
ple, a particle :need not impact, settle or diffuse the full distance re- 
quired by the equations, since its radius should have been subtracted. 
Here again this is appreciable only when d = 0.01. But then either 


the obtained fraction is negligible, so that correcting it is unneces- 
sary, or the correction becomes negligible when the fraction becomes 
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appreciable. At most, these corrections increase the per cent retained 
by only a few per cent. 


Calculations. We proceed to calculate the fractions of the air- 
borne particles which are retained in the various regions. Since we 
wish to compare results from deep breathing with shallow breathing, 
and in the former case the alveolar air exchange is a very large frac- 
tion of the tidal air whereas in the latter case the alveolar air ex- 
change is negligible, we shall proceed with a somewhat more detailed 
calculation than that carried out by W. Findeisen. From equations 
(1), (4) and (6), together with the numerical values in Table I we 
can estimate the fraction of particles removed, due to any one fac- 
tor, as the air passes from one region to the next. Then by successive 
calculations, which we shall now outline, the total amount removed 
in any one region can be obtained. 

From the values given in Table I, the velocity in each region 
and the time of passage are calculated. Since part of the time is 
spent before reaching the region in question, the corrected times, ra’, 

=a-1 
are obtained by multiplying 7.’ by the fraction passing, 1 aes Va. 

B=1 
These values are used throughout in the calculations. For the 4 sec. 
cycle, 300 cc./sec. the values are given in the last column of Table I. 

For any region a, the fraction P.* of the material entering it, 
which is retained in the region, is the combined probability due to 
the three main factors: Pat = Dg + Iat + Sat — Do'Ia* — Di‘iSa’ — 
TgiSqt + Da'Ia'Sai where each probability is a function of the transit 
time 7‘. Actually, at least one term is negligible and often all three are 
small, so that the calculation is not as involved as it appears. Let Pa° 
be a similar compound probability for expiration. These probabilities 
differ only in that the expression for J in the latter is Ia(Vau, Ra) 
instead of Ia(va., Ra) as in equation (1), because of the reversed 
direction of flow. Let P.* be the probability of removal of the par- 
ticle during the time that the breath is held constant, so that J—0 and 
the time ra in equations (4) and (6) is 7/8, T being the time for 
a cycle. Let fa be the fraction of the particles in a given volume of 
air which is still air-borne when it reaches region a. Then the frac- 
tion R,’ of the material breathed which is retained in region a dur- 


ing inspiration is 
12 
au fave: Dey Bs fa= (1—P,) ---- (1 — Pan), (7) 
B=a 


that is, a fraction of the particles must have reached region a (fa), 
must be removed (P.‘) and only the fraction passing 
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1:2 Oy 
( Di Vp V0) 
B=a B=1 


is involved. 


Similarly, for R.’, the amount retained during the period that 
the breath is held, we have 


R= fal — Pa’) Pa® Va= fasr Pa® Va, (8) 


that is, the material must reach the region (fa) but not be retained 
on entering (1 — Px‘), and must be retained during the pause (P.”) ; 
the relative amount involved is Va. 

The relative concentration Cg in region f after the pause is 


Cs,8 = fp (1 — Pe’) (1 — Po"). (9) 


After the volume of air has passed back to region § — 1, the concen- 
tration becomes fg(1 — Pg’) X (1 — Pg") (1 — P%s). If a volume of 
air reached a depth #, its concentration when returned to region 
y < 6 would be 


Cig — Cag) — Pepa) Cl fee ee (ee (10) 


The relative amount of the material which had reached a depth 6 but 
is retained in a on passing back out is the difference Cag — Cag. The 
total removed in a from all regions is then: 


£2 
Ra? = DY Vg(Ca+,g — Cag) « (11) 
B=a+1 

Since particles on the average pass halfway into a region in which 
they stop, we shall consider removal upon entering the last region, 
but not upon leaving it. In order to obtain the total retained in re- 
gion a we simply add RA,’ + R,” + R,°, and the total retention is 
R, => (R.' + Ro” + R.°). As a check we have the relation that the 
amount retained plus the amount coming out in the expired air must 
add up to unity. Thus it is useful to calculate separately the amount 
which is expired. The amount of air starting at region a (which had 
reached a depth a« at the end of inspiration) is Va and the fraction of 
particles remaining on expiration is C,.. Thus the fraction in the 
expired air is given by 
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1— Ryp= VaCr,a 


ad=1 


ave (Pt) (1 Pt) ee (1 PS) (1 =P) (12) 


l— P.8 ee. Cle Pe) = 


In addition to being useful as a check on the calculations, this ex- 
pression enables one to calculate the concentration in any part of the 
expired air, since one can obtain the sum over various intervals of 
a. The sum in any interval divided by the sum of the volumes for 
that interval gives the relative concentration. 

Special Assumptions. In calculating the various probabilities 
some special assumptions have been made, a few of which are some- 
what arbitrary. The velocity of the air entering the mouth at 
18L/min. was taken to be 100 cm/sec., while the velocity of air enter- 
ing the trachea (through the glottis) was in this case taken to be 150 
em/sec. Since the mouth is held nearly horizontally, a factor of 1.3 
was introduced into the argument of the exponent, e.g. cos y= 1. 
A similar factor was introduced for the pharynx. On the other hand, 
a factor of 1/3 is introduced for the trachea, which is nearly vertical 
in the experimental procedure. Since the alveolar ducts and sacs have 
very irregular surfaces which make settling and diffusion more ef- 
fective, a factor of 1.3 was introduced for H , H’ and IJ in the calcula- 
tions of D and S. In addition, because of the essential spherical na- 
ture of the alveolar sacs, a factor of 1.5 was introduced, since there is 
a greater specific surface and no angle y involved. Since the alveolar 
ducts and sacs expand with inspiration, slightly larger values were 
used when the tidal air was larger. The three regions were assumed 
to increase by the same fractional volume so that the total volumes 
were in proportion to the residual volume plus one-half of the tidal 
air. The effect of mixing the air which is exchanged with the ‘‘dead” 
air is not taken into account in the calculations. This would tend to 
increase the calculated values of the retention. 

The respiratory pattern chosen for the calculations is the fol- 
lowing: Inspiration, at constant rate, 3/8 of a cycle; pause 1/8 cycle 
(7*, no inspiration or expiration) ; expiration 3/8 cycle; pause 1/8 
cycle. 

Results. In Table II, the calculations for one set of conditions 
are given to show the steps which must be carried out to obtain the 
values sought. Columns 2 and 8 are obtained from the data of Table 
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I and equation (1), columns 4 and 5 from Table I and equation 
(4). Column 7 is the compound probability obtained from col- 
umns 2 and 4, while column 9 is the same as column 5 since J = 0 
during the pause and D = 0 for such a large particle size. Col- 
umn 8 is obtained from column 7 by equation (7). Columns 11 and 
12 are obtained from the preceding columns by the use of equations 


nN 
o : 
st Sslocooooocoonnte/0 
= meh nN 
|| 9 .& “ . 
EGCliagcloooseonrcocendto | 
oe oO rc N 
an es > 
tS o 
<a an IES OLMIS Bi! SS ar Ba 
oo |.8 
lanl 
o oO Q 
AS WOlLAnAMMIONODANMN KRW} 
oS ral crt mod lor) 
Oo 1D 
© ri 
S CSCloaonenoaonnocoooc od |a 
S AIH ANN 7) 
Be il 
3 
a = lon SSOSSSSAAN © @ & | 69 
mice = 
Si sgBleywoleoSoSooSoonnDyt%/a 
Bios] @ N i 
= n S| 5 
Eilasie ecoonn nt Otoanw 
n ~ 9 | 8 
FS] 83/4 Hmoonr to) 
: ra ri N 
25 De i aS& = 
Bailie oS 
z S DHMAHHADOSOSOSOS 
a Ss pf fen am oo tae = 
f& oD 
Ho 
ee o So [oo an) 
oz = NI SCOSSSSOSOSOANRN a 
ae |B 
eI ws | ey, O/OSOSSSCSCONNrETASO|A 
Be Oo 8 ak N ina) 
Baca 
Bllos|sgnjloooonnnoetaRol|s 
i) 2b /8 a oo 
es a iS 
BUS Sl\Ael/HHAN4oarHyornvMoly 
> ss Non a 
3° os 
HoiAiMrowowwocs - 
q S Re et & 
z 
) 
— 
5 o qyloCSCSCSSSOONNGHFS|S 
Thelie 
= oH SCSSSTSONNH SSO 
els Sly 2S a7 |S 
o-/3 
vo ie.) re 
+ a cooocoonttoRe old 
6 9 i-s 
S) 18 1) 16 | 
ge As SOTAWAP ATS 9 
3 ~ a 
SONOrevosd o |e” 
S oP om teed me me ed S 
oD 
re 3 
x 
S[ MN osteo OAS TAS 


54 HUMAN RESPIRATORY TRACT 


(7) and (8) respectively. Column 13 is the compound probability ob- 
tained from columns 3 and 4. Column 14 is obtained from column 6 
and from the last part of Table II (Cau,g — Cag) as indicated by 
equation (11). The table for Cag is obtained from columns 7, 8, 9 and 
13, using equations (9) and (10). The table for Cass,g — Cag is ob- 
tained from the table for Cag by subtractions. These values are then 
used for-column 14 according to equation (11). Column 15 is just 
the sum of elements of columns 11, 12 and 14. Column 16 is the same 
as the diagonal elements of Cag. The elements of column 17 are the 
products of the elements of columns 6 and 16. 

Table III gives the summarized results of the relative amounts 
of various sized particles deposited in different regions under selected 
conditions. It may be noted that for particles greater than one micron 


in diameter, the diameter may be replaced by \/p d, so that density 
is directly accounted for except for very small particles. For small 


TABLE IV 

RELATIVE CONCENTRATION IN VARIOUS FRACTIONS OF THE EXPIRED AIR 

Fraction of Particle Diameter in Microns 

expired air 20 6 2 0.6 0.2 

300 cc/sec., 4 sec. cycle, 450 cc tidal air 
0-1/8 50 95 100 100 100 
1/8-1/4 20 50 89 98 99 
1/4-1/2 0 8 62 90 91 
1/2-3/4 0 0 30 70 70 
8/4-1 0 0 2 66 65 
800 cc/sec., 8 sec. cycle, 900 ce tidal air 
0-1/8 25 65 92 99 99 
1/8-1/4 0 4 44 85 82 
1/4-1/2 0 0 8 64 65 
1/2-3/4 0 0 1 sy 60 
3/4-1 0 0 1 57 60 
800 cc/sec., 12 sec. cycle, 1350 ec tidal air 
0-1/8 10 50 80 97 97 
1/8-1/4 0 0 20 75 70 
1/4-1/2 0 0 5 55 56 
1/2-3/4 0 0 5 54 54 
3/4-1 0 0 0 33 36 
1000 cc/sec., 4 sec. cycle, 1500 cc tidal air 

0-1/8 15 40 84 98 98 
1/8-1/4 0 0 57 82 82 
1/4-1/2 0 0 oe 75 15 
1/2-3/4 0 0 20 70 (ol 
8/4-1 0 0 6 56 55 
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particles of about 0.2 micron, the effect of settling is comparable to 
that of diffusion. The introduction of d’ = \/p d as a new variable 


introduces d into Sa, but V1/d =V/ V/p/d’ into D,. Thus even for 
p = 2, p+ = 1.2; the error when both terms are comparable is about 
10%. If the calculated retention is 30%, the correction would be only 
about 3% for d = 0.2 microns. 

Table IV gives the relative concentrations of material in the vari- 
ous fractions of the expired air. These values are calculated from the 
values of Va Cia. For example, the first entry in Table IV is ob- 
tained as follows: In Table II the sum of Va Cic for all a in the first 
1/8 of the expired air [Mouth + Phar. + (0.045/0.06) Trachea] is 
0.022 '+ 0.021 + 0.022. Thus the concentration is .065/(1/8) = 50%. 

We next compare the results of these calculations with some 
available data (cf. Landahl and Herrmann, 1948). For corn oil drop- 
lets, the calculated values are somewhat too high in the neighborhood 
of fifty per cent retention, whereas for methylene blue particles, the 
calculated values are somewhat too low in the same range. For cal- 
cium phosphate dust, taken to be particles composed of packed 
spheres, the calculated values are appreciably higher than the ob- 
served values. 

If sodium carbonate remained dry, its density should change the 
effective diameter by a factor of 1.5. But if the particles were all at 
equilibrium with completely saturated lung air, they should take up 
water to a diameter of about five times the original. Since the large 
particles may not have reached equilibrium, the effective factor to be 
applied to the diameter may not be constant. If a factor of two and 
one-half is applied, there is fairly satisfactory agreement at 18 and 
60 liters per minute for the four second cycle, but not for the other 
cases. The comparison in these last cases is perhaps unjustified be- 
cause of the hygroscopic property of the substance. Therefore, for 
the present, we may conclude that there is no serious discrepancy be- 
tween these theoretical results and experimental data. 
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_ Some thermodynamic aspects of steady systems are considered. The 
time rates of change, “flux,” of various thermodynamic quantities are 
formulated. In particular the free energy flux in the steady state, the 
difference between the free energy flux in the steady and time depend- 
ent states and the change in free energy flux upon transition between 
steady states are discussed. Equations are derived which exhibit the for- 
mal similarities and differences between the free energy flux and the 
conventional free energy change. The temperature dependence of the 
steady state rate is examined and conditions for “mastery” by a single 
step discussed. A brief discussion of the role of rate in the coupling of 
exergonic and endergonic reactions is given. 


Introduction. In two previous papers (Hearon 1949a & b, here- 
inafter referred to as I and II respectively) some aspects of the ki- 
netics of steady state systems were discussed. In this work two mat- 
ters, both of importance and interest, were not considered viz., the 
time course of the attainment of the steady state and the energetics 
of steady state processes. The former will be reserved for a subse- 
quent paper now in preparation. The latter, which will be considered 
here, falls not strictly in the realm of formal kinetics. As will be seen, 
the problem of the energetics of steady state processes is intimately 
connected with the rate of the process. Thus the separation of the 
analysis into kinetics and energetics, while demanded in general by 
the strict adherence to the mechanical analogy, is logically faulty 
and of no use in practice. 

We shall be interested here not only in computing the usual ther- 
modynamic quantities for systems in which processes are proceed- 
ing at a finite and specified rate but also in the irreversible or dissi- 
pative aspects of such processes. The reasons for this emphasis are 
many. The computation, for example, of the irreversible production 
of entropy is of importance in principle and practice alike as has been 
discussed by R. C. Tolman and P. C. Fine (1948). The concept of 
irreversible entropy production or of the “uncompensated heat’ of 
Clausius (De Donder, 1936; Van Rysselberghe, 1935; La Mer et al, 
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1949) permits the second law of thermodynamics to be expressed as 
an equality. This is clearly an extension of the conventional formu- 
lation of some importance. The inclusion of the irreversible entropy 
in the first and second laws of thermodynamics is also of great prac- 
tical importance from the standpoint of the discussion of the efficiency 
of thermodynamic processes as has been discussed in some detail 
(Tolman and Fine, 1948). ; 

With the knowledge of the irreversible entropy production one 
can calculate in a simple manner the total rate of energy dissipation 
(a concept which is made precise in what follows) in a system in 
which irreversible processes are occurring. This quantity, in quali- 
tative disguise, has entered biological theory at an early date. 

Long ago A. J. Lotka suggested that “ ... . natural selection 
tends to make the energy flux through the system a maximum, so far 
as compatible with the constraints to which the system is subject.” 
(Lotka, 1922.) This postulate and its logical basis are also carefully 
stated, respective of the restrictions that the. supply of energy and 
matter not be limiting, in Lotka’s book (Lotka, 1925) and have been 
discussed more recently elsewhere (Lotka, 1945, 1948). The term 
energy flux is defined explicitly (Lotka, 1922, p. 150, footnote 5) as 
“... the available energy absorbed by and dissipated within the sys- 
tem per unit of time.” (this author’s italics). 

N. Rashevsky (1943a) gave analytic form to the maximum prin- 
ciple of Lotka as a variational principle and subsequently (1943b) 
introduced a general Hamiltonian principle designed to care for the 
events of classical physics as well as those of biology. In application 
of this principle to the theory of amoeboid movement (1943c), cell 
division (1943d) and organization (1943d) by Rashevsky the term 
energy flux has taken on the meaning of energy exchange between 
biological units or between biological units and their environment. 
The quantity to be maximized is the total positive flow of energy. 
Dissipation of energy due to chemical reactions* is considered as a 
negative flow; flow from the system. Dissipation due to viscous forces 
which depend upon (generalized) velocities is cared for separately 
in the Hamiltonian principle (Rashevsky, 1943b). This would seem 
to be a rather different interpretation of flux of energy than Lotka’s 
original definition. Furthermore in his most recent discussion Lotka 
(1948) introduces a slight change in point of view stating that 
“,... those species have the advantage which know best how to take 
and put to favorable use the stream of energy available in nature. 

*The dissipation term used by Rashevsky to account for ( spontaneous) chemi- 


cal reactions is quite different formally and conceptually from that discussed in 
IV following (p. 86). 
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Collectively this will tend to increase the flux of energy through the 
system, which flux may thus be said to be approaching a maximum.” 
(Lotka’s italics). This implies that efficiency plays a definite role.* 
The point of this somewhat lengthy discussion is to exhibit the 
importance, in general biological theory, of the dissipative aspects of 
energetics. The question of the relative roles of energy flow or ex- 
change of dissipation and of efficiency in an adequate maximal prin- 
ciple, if such be possible, is an involved one: we eschew it here. It is 
pointed out, however, that the question of the possibility of such a 
maximal principle in biology was a factor in the initiation of the 
present investigation. Further we submit that the thermodynamically 
acceptable form for the quantities required for such a formulation, 
l.e., dissipation, rate of change of free energy, etc., are available in 
the results of this paper and the more general formulations of irre- 
versible processes to be referred to in IV. 

The results to be considered divide naturally into two groups: 
1) The formulation of the various energetic quantities for systems 
which are homogeneous and in which chemical processes are occurring 
at a specified rate. These results are presented here. 2) Similar for- 
mulations for non-homogeneous systems (systems in which there are 
spacial gradients of composition), for geometrically confined regions, 
such as the cell, and their associated diffusion fields. These are pre- 
sented in part IV. The above corresponds to the division of purely 
kinetic topics into those discussed in I and II respectively. 


Closed, Homogeneous Systems. 


In this part we will discuss the general set of consecutive reac- 
tions considered in I and II; familiarity with these papers is pre- 
supposed. We regard the system as closed and hence exclude 
the interchange of matter between the system and environment. 
The constancy of the concentrations of the reactants, [C;], and of 
the products, [B;], is thus supposed to be cared for in one of 
the ways previously discussed. The manner in which the follow- 
ing problems are formulated in general will depend upon the in- 
formation available and upon the conditions imposed. For any given 
set of conditions it is well known that there exists a characteristic 
function giving rise to what J. W. Gibbs (1928) termed a fundamental 
equation. It is also well known that knowledge of this characteristic 
thermodynamic function and the fundamental equation allows the 
derivation of all other thermodynamic properties of the system. The 

*That Lotka had not previously ignored the factor of efficiency is evident in 


an addendum to the 1922 article. The reader is referred to this and to Lotka 
(1925) for his arguments on this point. 
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conditions of most general interest are those of constant tempera- 
ture and pressure. The point of departure will thus be to choose the 
Gibbs’ function or free energy, F', as the characteristic function. Fur- 
thermore for isothermal processes the irreversible entropy produc- 
tion is most simply obtained from the free energy change. It will be 
necessary to assume as available extra-thermodynamic information 
in the-form of the rate constants and their temperature dependence. 
Otherwise, however, the formulations are entirely within the frame- 
work of classical thermodynamics, a feature which is discussed at 
the close of this paper. 

The free energy will be taken as a function of the temperature, 
T , the pressure, P, and the composition which is specified by 7; , 


7=1,2...., the number of moles of the 7th constituent. The fa- 
miliar differential, at constant T and P, is 
dk => Und, , (1) 
k 


ena Aaa 0: | 
i} 


aa 


j #k. They are 


T,P,n; 


oF 
where the mu, are the chemical potentials( ued ) 
Ni; 


taken to be of the form 
KK ua + RT ln Ci 4 (2) 


where u,° is the standard chemical potential of the kth solute, a func- 
tion only of JT and P when the solvent is specified; C; is the concen- 
tration of the kth solute. Relation (2) embodies the usual assump- 
tions (Guggenheim, 1930, 1936) regarding perfect dilute solutions; 
in what follows the modifications necesary when these assumptions 
no longer hold will be noted at the appropriate time. 


Time Rate of Change of F. 
For the chemical reaction 


— »,R, — mR, =3P3 + %4Ps (3) 
one may, following Th. De Donder (1936), define the quantity & by 
dn; = yi aé ‘ (4) 


The »; are the number of moles of the ith constituent produced in the 
reaction (3) as written. The quantity & is the extent of the reaction 
(see also Mayer and Mayer, 1940; MacDougall, 1939). If the dn; 
from (4) are substituted into (1) we obtain 

oF =y 

ae = : Vi Uke (5) 


The time rate of change of the free energy is given by 
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a =(2)%- dé 
in ae Gal l"a) S- (6) 


ee 
The quantity He is the rate of the reaction in the sense that 


d& d(ni/v:) 
dt dt 


Thus dé/dt is the time rate of change of the number of equivalents, 
ni/¥i of any constituent and is easily expressed in terms of more 
familiar chemical kinetic expressions. The quantity S| », um, is the free 


5 (7) 


k 

energy of the reaction. It is frequently written simply as AF to desig- 
nate the change in free energy, per equivalent, of the reaction; it is 
actually the slope at any point of the F vs. & curve. When the ;, from 
(2) are substituted into (6) the concentration values must then be the 
instantaneous values. If several reactions are to be considered, a é; is 
defined for each reaction and the equations corresponding to (5) and 
(6) are obtained by summing over 7. Implicit in (5) and (6) is the 
assumption that the mechanism of the overall reaction is correctly 
portrayed by (38): There can be no intermediate steps or reactions 
and the disappearance of y, and v. moles of R, and R, must be strictly 
accompanied by the appearance of v; and », moles of P; and P.. 

The expression such as (6) for the set of reactions discussed in 
I may be derived, taking into account that there are many steps and 
that every X; participates in at least two such steps, as follows: 
Since every |v;| is one, there is no distinction between equivalents and 
moles. Consider the ith step 


Xin +O, 5 X%;,4+ B;. 


The free energy change per mole is AF’; and the rate (moles/volume/ 
time) is v;. The rate of change of the free energy for the ith step, 
per unit volume, is thus AF; v; and for the complete set of reactions 


dF 
—=F(t)=VSARY, (8) 
at k 


where V is the total volume of the system. The dimensions of F(é) 
are energy /time and the units will usually be calories/sec. Denote by 
[Xz] the concentration of X; when the system is in a time-dependent 
state and by [.X;] the concentration when the system is in the steady 
state. (Note that in I [Xx] signified the steady state value with no 
qualifying notation for only steady state systems were discussed 
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there.) Using the notation of I with the above modification, 
w+LXx] 


AF, = RT Ge tat (9) 
where by definition, 
[Xo] = [Xn] =1. (10) 
For v; we have 
Vy = WeLXnal — wW-«lXx] . (11) 


Expression (9) is easily seen to be identical with [recall equations 
CEy ors) | 


AF, — Be. = (eee (uc, Br Mx, ) 
By [Xx] ee 
— ARS 40RT Ine 
[Ci] [Xi] 
when it is recalled that 
Me? r Mx, i, (uo x i ux, °) os AF’,° | 
(13) 


——RTiInK;,, | 


and that the w’s are given by (2). Equation (8) with AF; from (9) 
and v;, from (11) yields 


w«e[ Xx] 
Wr[Xxa] : 


Equation (14) correctly accounts.for the fact that every X; is pro- 
duced in the kth step at the rate v; and is consumed, at a different 
rate, Ve, in the subsequent, (k + 1)st, step and that although the 
production and consumption occur at the same instantaneous value 


of [X;] the AF’s of the kth and (k + 1)st steps are in general differ- 
ent. 


F(t) =VX (wil[Xea] — wa [Xa]) RT In (14) 


Rate of Change of F in Steady State. 


When the system is in the steady state the value of F is no longer 
time dependent and is denoted as F,. Further for every k 


[Xx] = La]; 
(15) 


V,— V3. 


Equations (11) with (15) lead to the set (2) of I and 2, is the steady 
state velocity discussed in I and designated there as v. 
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Equation (14) then becomes 


F.=V(SRrIn— Jo, 


k=1 Wk (16) 
el 04 I ETP 
Equation (16) is readily verified by noting that 
n [Xx] n n 
Sin =m| (xa - | 
wig. oee pts SAL Nes] 
n n (17) 
=mn| 1X) 7 |=o, 
ket jor LX] 


where in (17) the definitions (10) have been used in effecting the 
change of index. It is clear that in (16), AF, the free energy change 
of the overall reaction, is, because of equations (1) of I and (17): 


AF =>" AF, 
k 

(18) 

Bul tB. eee Be 

= AF + RT In ee Bal [Bol 1182] Le 
[Co] [Cy] Aion [C,] 
where 
AFPo=—RT In K=—-RT In ILE, 
k= 

(19) 


=> us’ —>d Me,” 
k - k Md 


which is the conventional but less useful form. So long as the system 
is in the steady state the change in free energy per unit time has the 
constant value F,; it can be shown in a simple manner that this is 
always a negative quantity. Equation (16) is written, for this pur- 
pose, in the form 


II w+ ITlw,— Il Wk 

pav( arms > ees ). (20) 
I wy a 

k 


The last factor in (20) is v, from equation (3) of I with ¢ as de- 
fined in equation (18) of II: ¢ is a sum of positive terms [eq. (21) 
of II]. The second factor in (20) is AF from (16). If 

Iw, > wz ,v;>O0and AF <0, ~ (21) 


k k 
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and if 


Tw, < WZwsz,V’s<Oand AF > 0, (22) 


k k 


thus F, < 0 under all conditions, the equality holding for the equi- 
librium condition. Similarly F(t) is always negative, for by the same 
reasoning every AF; v; in (8) is negative as can be seen from (9) 
and(11). In words, if the first inequality of (21) is satisfied, the 
steady state rate is positive and the last equality shows that the pro- 
cess is spontaneous. If the first inequality of (22) is satisfied, the B; 
are converted to the C;: AF for the reaction as written is positive, 
but the process proceeds in the direction of spontaneity for the re- 
verse of the reaction as written takes place. From the above it is 
clear that AF = 0 is a sufficient condition for v,; —0. For the sys- 
tem discussed it is also a necessary condition; the vanishing of vs; 
necessarily implies that AF —0. In general however this last state- 
ment is not true. Is should be made clear that although we have re- 
garded every [C;] and [B;] and hence every w; and w-; as constant, 
the conclusions F, < 0, F(t) < 0 drawn from (21) and (22) do not 
depend upon this condition. What we have done is to admit the net 
conversion of the C; to the B; but not admit changes in these concen- 
trations. If one considers an open system by requiring that the [C;] 
and [B;] are maintained constant in the system proper by flow from 
and to the environment (see discussion in II), then in the steady state 
there is no net production of any C;, B;, or X;. For such a system 
the condition, cf. the second equality of (18), that the difference in 
chemical potentials of the products and reactants (which we will call 
the potential difference) in the system is zero* is a sufficient condi- 
o[Cs] sel Bile kal sae 
SS ee a ee ciel 
ot ot ot 
necessary condition for every rate may vanish without necessarily 
implying that the potential difference is zero. For X;, 
dL x 
= . = (wel Xen] — we LXa]) — (Wee [Xa] 


tion for the vanishing of every rate, 


(23) 


— W- (K+) [Xi] ) = Uy — Views; 


the rate vanishes if v; and v;., vanish separately: This is the condition 
of equilibrium and the potential difference is zero. The rate also van- 


_*Note that this condition can be attained when and only when the (fixed) 
environmental concentrations are such that the equilibrium condition may be sat- 
isfied. Similarly the closed system discussed here can come to equilibrium only 
if the w’s are arranged to satisfy the equilibrium ratio unless the condition of 
constancy be relaxed for some of them. : 
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ishes if Vi = Vir, this is the steady state and the potential difference 
is not zero. 

It is in fact true that when such an open system is in the steady 
state every extensive property of the system is time independent. 
There is no time rate of change of F in the system but only in the 
environment. This is of considerably more importance in the topics 
discussed in IV. 

C. Eckart (1940, see addendum p. 924) has given a rather de- 
tailed proof of what is equivalent to our propositions F (i) =. 0s 
F, < 0 and discussed whether the usual kinetic rate formulae are in 
accord with the general thermodynamic principles discussed in his 
paper: Both the proof mentioned and demonstration that the formal 
kinetic rate laws are in accord with general thermodynamic princi- 
ples are based on expressing any concentration as 


[Ci] =exp { (ui — ui°)/RT} . (24) 


Inherent in the rate laws used by him is the fact that 
> w°=—RTInK; (25) 


cf. equations (13) and (19) of this paper. It seems necessary to only 
mention this here especially in view of the fact that the absolute re- 
action rate theory (Glasstone et al, 1941), as opposed to classical col- 
lision theory, brings kinetic rate theory into complete accord with gen- 
eral thermodynamic principles. Later it will be pointed out that De 
Donder (1936) has already given a theorem which will be shown to 
be equivalent to F < 0. 


It is worth mentioning that the set of reactions discussed in I 
and II is not very general in its stoichiometry, i.e., every |v;| = 1 as 
mentioned. The modifications of what has been formulated here for 
a different stoichiometry are simple and obvious. The reactions can 
be made as general as is necessary: There may be any number of 
B’s or C’s in any step or there may be two of a given C;, in which 
case the square of [C;] appears in w;. The only restriction is that 
each step be first order with respect to the X’s involved, for only then 
is the steady state solution given in I valid. Be this as it may the set 
of reactions was intentionally set up in this manner. For it is re- 
garded as most improbable that a single step be of more than second 
order. It is well known that three-body complexes are formed in many 
catalytic systems of biological importance but they are formed step- 
wise by first and second order reactions. 
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The Difference F(t) — F,. 


It is of some interest to inquire as to possible extremal proper- 
ties of F,, i.e., whether F under steady state conditions is always 
greater or always less than under variable state conditions. To an- 
ticipate a later discussion it may be stated that this will also answer 
the same question regarding the rate of dissipation. It can be shown 
from (14) that the 0F/é[X;] vanish simultaneously only if the con- 


ditions of equilibrium, v, = 0 for all k, obtain. The difference, 
F(t) — F, which may be positive, negative or zero is given by 
bt AP eile (26) 
k 


Equation (26) results from (8) and (16); the AF; are to be taken 
from (9). The result may be expressed as in (26) because the addi- 
tion, to (16), of RT times (17) changes nothing whether or not the 
[X;] are steady state values. Certain forms of (26) are of interest. 
Substitution of (9) into (26) gives 


F(t)—F, °% Wx aes [X%] 
1 =3( Brin ) co <0) +3 (RT In ) me. on 
V ke1 Wk : k=1 [Xia] = ! : 
Because of (10) the last summation in (27) can be written 
> (RT In[ Xi] ) (V4 — Ver) - (28) 
k=1 
Now the rate of change of [X;] is, from (23), 
d[ Xx] 
ve = We LXra] — (Wa + Wier) [XE] 
29 
ae W_(k+1) [Xia] ’ ) 
or, in terms of the displacements from the steady state values, 
dAX;, 
ee = Wy, AXi-a — (wa + Wri) AX, 
Pie W_(k+1) AX aes = AXs ’ gh 
where Lt | Se aa 
AX, = [Xx] — [X'] . (31) 


But from (23), (29) or (30) is an expression for ye Nee Also. 


Un — Vs = Wr AXk-a SW AX; 2 (32) 
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Putting equations (28), (30), (31) and (32) into (27) gives 


F(t) —F, n ey 
sista = 2 ( RT in =| (w. AX 54 uy bias AX) 
V Heal Wk (33) 
+3 [RT In( [Xi] + AXx)] AX:. 
kei 


The definition (10) gives, of course, 
AXG=AX> =0: (34) 


Equation (33) involves only steady state values of the [X;] obtain- 
able from equation (6) or (7) of I and the AX;, obtainable from the 
solution of the set (30). Because of (34), explicit introduction of 


AX;, from (30) into (33) gives* 


Foy—f. -* 
SAA (WAP + Wie AP pi) (35) 
V k=1 
in which AF;, is to be taken from (9) and, of course, 


AF,=— AF; . (36) 


By (81), (35) can be expressed in terms of the [-X;°] and the 4X;; 
as it stands it involves the time dependent X-concentrations and the 
AX,. As a simple example the case A S X, S B will be discussed. 
Using the form (35) 


F(t) —F, 
V 
[X,] 


AX, (er are )+ is Ax,( RT a aoe ) (37) 


1 
The solution of (30) for 4X, is 

AX a Xe (38) 
where a =k. + & and AX,° = [X,]° — [X,’] is the initial displace- 


ment from the steady state. Equation (37) vanishes for very large 
t when AX, — 0; it also vanishes for 
Ke {A} 2/4 
Exe] =e al ieee (39) 
(K,/[B])"/” 
From (39) it is seen that 


*Note that (35) may be formed by taking the “self term” of (29), multiply- 
ing each w by the corresponding 4F' and summing. 
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if Kix Pe Ke 9 [X,] lgrd Ky [A] ; 
if ka<<hke, [Xs] ~ [B]/K. 


Under these conditions (37) vanishes when the first or second step 
passes through the equilibrium position. For the [X,] values in (40) 
are approximately equilibrium values but clearly the two conditions 
in (40)-are mutually exclusive. The time at which (387) vanishes, 
other than t = o ,is 


(40) 


n : (41) 

Gi tANae 

From (38) 4X, and AX,° must have the same sign, further |4X,°| “ 

|AX,| , thus ¢ from (41) is real and finite, for 1/a is inherently posi- 

tive, and at this time (37) vanishes as well as at very larget. 

For more complex cases the 4X; from (30) are of the form 
AX; => kj et, (42) 
4 

and the analysis is not so simple for the 4X; themselves need not be 

monotone in t (Burton, 1939) a matter which will be discussed in V 

of this series. It is of interest to note for the simple system discussed 

here that “overshoot” in a thermodynamic quantity is possible when 

the concentration [X,] is capable of no “overshoot.” 


Change in F,, Upon Transition 


The expression for the change in F’, between one steady state and 
another will depend upon what has induced the transition between 
steady states. If any of the [B;] or [C;] be changed both v, and the 
potential factor, what we have called the potential difference, given 
by (18), will usually change; only one part of the potential factor 
will change for AF° is independent of concentration. If one imagines 
one or more rate constants to be changed, v, is changed and only AF° 
of the potential factor is changed. Finally if a given C; is conserved 
in the system, such as is true of catalysts (see 1), then there is some 
Bn,m > 7, which is identical-with C;. The concentration of this com- 
ponent may enter into v, in a complicated manner ( cf. I) but its con- 
centration always cancels from the potential factor as is obvious from 
(18). Thus if the concentration of a catalyst be changed, only v, is 
changed. These cases may be indicated as follows: In (18) denote 
the concentration ratio in AF by R,. Thus R, is the “product to re- 
actant ratio” for an n-step system. Then (16) becomes 


F,=V[AF°+ RT mu R,] 2, . (43) 
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If any non-catalytic [Bi] or [C;] be changed: 


F > nome 1 dere 
aoe = AF? (Vg. — V5) + Vs2 RT In Rus 
(44) 
=O) tul or teas « 
If any rate constant k; or k_; be changed: 
Fy Tay Fs, 
7 pas AF (2)° Ve — AF n° ta + Oe 
(45) 
ee) tel tbe we 
If any catalytic or conserved [C;] be changed: 
P50 a FS 
pS AP! + RT ltt Rn) (V2 — Ver) « (46) 


Equations (44), (45) and (46) give the excess of rate of expenditure 
of free energy, per unit volume, for the second steady state over that 
for the first steady state when the transition is induced by the agency 
indicated. The subscripts 1 and 2 refer to the respective steady states; 
quantities with no such subscript are identical in the two states. 
Equation (46) applies not only to the change in concentration of 
enzymes but also to the change in concentration of many metabolites 
which function cyclicly and are conserved. The role of such cyclic pro- 
cesses has been discussed in great detail by H. A. Krebs (1948). The 
difference in free energy between the two states may be obtained by 


tz 
F.—Fi= { F(t)dt, (47) 
ta 


where t, is the time at which the transition was initiated and ¢, is a 
time chosen such that the new steady state is approached with any 
degree of approximation, since in principle the transition takes for- 
ever. It seems of doubtful value to attempt to carry (47) further in 
any general case. For that matter it is the heat change between two 
such states which is experimentally amenable and of interest. 


Time rate of Change of Hand S. 


Thus far we have confined the discussion to the quantity F. A 
thorough treatment of the subject, some aspects of which are being 
discussed here, would exhibit the derivation for a steady state system 
and for a system approaching such a state the counterpart of all the 
usual qualities discussed for equilibrium reactions or for other pro- 
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e 
cesses, not withstanding the time factor. Actually not all counter- 
parts exist; of those which do, not all are of value or interest. As 
mentioned it is a familiar procedure to derive from a characteristic 
function and its fundamental equation the remaining thermodynamic 
functions. 

It will be shown here what results when this same procedure is 
adapted formally for the function F,. This will! lead naturally to the 
time rate of change of other thermodynamic functions and will raise 
in a definite manner the question of the temperature dependence of 
v,. This latter question will be treated. In the first place it is evident 
from (16) that if 

aLf (a) -F] _ 
0x 


Y being the partial derivative of the product of any function f(x), of 
some parameter x, and the free energy F with respect to the para- 
meter x , then it is true that 
O[f(%) - F] 
ox 


yo (48) 


olnv; 
—=Y+ f(x) -F——. (49) 
Ox 


The quantity Y is 


gas ViAAcy, (50) 
a Us, 

dt 

where AY is the change in Y for the overall reaction; Y will here 
be some extensive property of the system but this need not be true in 
general. For want of a better terminology we will call Y the Y-flux. 
Thus F, H, S, E, etc. are the free energy-flux, enthalpy-flux, total 
entropy-flux, energy-flux, etc. Some examples of (49) are 


oF./T HH F,dlnv, 


II 


Y 


aT pet pe apes (By) 
oF, ETS ys poet Ue, (52) 
or 
oe y = FS In Vz ia 
oP oP 


The quantity AV is the change in volume accompanying the overall 
reaction 


AV=3.Vs,—3Ve, - 


q=1 i=0 
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In taking V the volume of the system as fixed we have regarded AV 
as negligibly small which it, in fact, usually is though it is determin- 
able from the P-dependence of F. An alternative form of (51) is 


Of s/ Les F,olnv, 
01/T Teo Tas 


the derivative in the last term of (54) is a familiar one in the study 
of rates in their temperature dependence. 


The total entropy-flux and enthalpy-flux may be discussed briefly 
as follows. By the general definition (50) 


(54) 


S=VAS2,. (55) 
By carrying out the differentiation indicated in (52) it is seen that 
S=V(AS°—RinR,)2;, (56) 
where RF, has the significance discussed in connection with (43) and 
Bg Sear. (57) 
fe 8 


By using the familiar 4F° = AH® — TAS® and (19), equation (56) 
may be put in the form 


A H° Wk 
s=v(——-sRn=)»,, (58) 
T k Wk 


comparable to (16). Similarly 
H=VAHdHv, (59) 


gives the heat absorbed per unit time. The heat absorbed in the tran- 
sition from one steady state to another is given by 


te ¢ 
Fie ae { H(t)dt, (60) 
ti 


where the limits of the integral (60) are as discussed for (47). It 
would be of great interest to investigate the usefulness of (60) in 
problems such as “heat of shortening” (Hill, 1949) in which the 
situation is exactly the transition from one steady rate of heat evolu- 
tion to another. Because we have considered perfect dilute solutions, 
the values of AH and H are identical with the corresponding stand- 
ard quantities A4H° and H°, but we have not used this fact. If activ- 
ity coefficients are inserted into the w’s the following change must 
be made: there will appear in (56) the additive term 
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oT 


where y, is the ratio of activity coefficients formed as in Rn»; in (58) 
AH will appear in place of 4H°; AH in (59) will be 


(61) 


? 


0 ln Yn 
AH*-+- R : (62) 
aphagk 
There is one further derivative which should be mentioned, viz., 
7 Inv; 
Salat ae J Te Ty a ness (63) 
o[Ci] TY 
In general with v, from (3) of I, (63) is not simple. From (16) it 
follows that if any wz, —=—0, AF =— o. It is possible however for 


the overall reaction to be irreversible so far as vs; is concerned and 
still have AF’ maintain its form. In short II w;, may be neglected 


with respect to in w+ and AF by no means bé infinite. In fact the 


reaction can be eoteider ed as irreversible for values of the w’s which 
permit no terms to be dropped from AF. For example then let 
W_, =0. Taking v, from (24) of II there is obtained 


ale, = ®; VU; 
@LCi] Vn‘? [Ci]? 


Eliminating 6; , because it is less simple to construct (see II) than 
Vi. by using (24) of II, (68) becomes, 


@ Eiaeae| 
eLci] = [Ci] 


(64) 


of Vn? VO, AT) (65) 


T-dependence of v;. 


In (51) and (52) and (54) there arises the temperature deriva- 
tive of v,. Having given (56), (58) and (59) for S and H it remains 
to complete (51), (52) and (54) by exhibiting explicitly the T-de- 
pendence of v,. This is in itself an important problem in steady state 
kinetics and it may be anticipated that the expression will involve a 
large number of thermodynamic quantities. It is commonplace in 
biology to assume that the T-dependence of v, is determined by the 
thermal characteristics of the slowest step in the sequence provided 
only that it is sufficiently slower than all other steps. This assump- 
tion is tantamount to the assumption that the expression for the over- 
all rate, v, , is that for the slowest step. This assumption is true only 
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if other conditions are met as is apparent from the discussion in I; 
further, it is by no means confined to biology (Glasstone et al, 1941, 
199; Branch and Calvin, 1941, 375). It will be seen that in general 
the T-dependence of v,; involves the standard heat of the overall reac- 
tion, 4H°, the standard heats, AH;,°, of the individual steps and every 
standard heat of activation, 4H;,i, of the forward rates. The general 
expression is unwieldy and it is a separate problem, pursued to no 
great extent here, to assess the various permissible approximations. 
We begin by dividing numerator and denominator of equation (3) of 
I by IT w,, obtaining v, in the form (66) which contains the maxi- 
k 


mum number of equilibrium constants and the minimum number of 
individual rate constants all of which refer to forward steps, viz., 


1— R,/K 
Cae eS (66) 
> Ax 
k=1 
where R, and K have already been defined, [see (18) and (48)-(46) ]. 
The A; are additive terms, T;,, [equation (21) of II] of ¢, equations 
(18) of II and (20) of this paper, divided by IT w;,. They may be 
k 


written as 


A =(mr-)( 72 )- 
: j= jor LCs] J We 


Rin 1. (67) 


Kanes Ws : 
Kura = 1 if ii 1 ’ 


where R,-, is the “product to reactant ratio,” constructed as R,,, for 
the first k — 1 steps and K,,., is the overall equilibrium constant for 
the steps 1 through k — 1 as is clear from (67). From (66) are ob- 
tained 
0 0 
(s Ax) — (—R./K)+ ( R./K ae 
k ri Oak 
ee bee -—$_——————_,, (68) 


fy DA Ag 
jk 


fo) 


and 
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0 0 
omv, oT ini) OL : 


= —_______ — —_—__, (69) 
pul tg > Ax 
Kk 
Equation (68) or (69) may be simplified as follows: 
fy fern” 
oT oe FOU fey bes 


where AH® has been defined and we have used the well known rela- 
tion, 


Oln K 
7 =A H/RT. (71) 
Also 
OAy Ria 0 1/Ka 7-1 may Rix alee ln Wk ; (72) 
0 T Wr 0 T Kans Wr 0 it 
and from (71) 
0 1 K, =] A He ei 
t g = = ? (73) 
0 T Kana RT? 


where A H°,;, is the standard heat of reaction for the steps 1 through 
(k — 1) and in view of the make up of K,-, from (67) is, of course, 


k-1 

AH x.=SAHS, k>1 
j= 

=) igeale 


Finally since w, is a forward rate constant multiplied by concentra- 
tion terms, with the usual approximation,* (Glasstone et al, 1941) 


(74) 


0 ln Wy 


=AH,/RT, (75) 


where AH,‘ is the standard heat of activation of the forward kth 
step. Using (73) and (75) in (72) gives 
0 Ax | AH y.-+ A Hie 
— = — A, ) ———__—__ } . 


oT | RT? .75) 


kT . 
*Strictly k; = ( br ) exp { — AF ,t/RT} with k the Boltzman and h Planck’s 


constant, the variation of In T relative to 1/T is neglected. This is usually done 
(Glasstone et al, 1941). 
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Substitution of (70) and (76) into (69) gives 


dlnv, —AH°/R 2 Ax(A 1,4. + A Ht) /R 


ly ae Fetes ey sa 
k 


(77) 


It is easily seen that if the overall reaction is irreversible, say, be- 
cause w_, ~ 0 then (66) gives 


Ve=1/> Ax, (78) 
k 


and (77) becomes 


; — 3 An (4 Wn. + A Ait) /R 
Olt Vspagee x 


a1/7 * Sid, (73) 
k 


There are several points to be noted. If the overall irreversibility re- 
sults from any w; being large then all A, with k > 7 vanish from 
(79) and (77), if it results from any w_; being small then all A; 
with k = 7 + 1 vanish. If a given w; is small then all A; with k < j 
will be negligible. These are conditions on the w’s and not on the 
rate constants; they may be met by the proper magnitude of the rate 
constant or concentration term involved. If a given k; is small then 
the A; with k => 7 will predominate; the value of A Hi? may be large 
relative to all others but the products A, A H;* with k > 7 are cer- 
tainly not negligible relative to A; A H;?. Further it is not always 
lawful to assume that a large k; implies a small A H;# or the converse 
for there are known cases of relatively rapid rates associated with 
strikingly large A H? (Glasstone et al, 1941); the apparent discrep- 
ancy being accounted for by the entropy of activation AS? which oc- 
curs in AF? (cf. footnote p. 74). Finally it must be remembered 
that the ‘A H°, ;., from (74) are not independent of the A H;+ since 


AH,=AHi—A Hz. (80) 


It is clear from the above and from the discussion in I that there 
are many ways of achieving overall irreversibility and conditions of 
‘limiting rates’: 1) If some w; is large, then as discussed in I, if 
j #1, the expression v, reduces to that for the first 7 — 1 steps. 2) 
If some w-; is small, then the expression v, reduces to that for the 
first 7 steps. These statements are general so long as the w in ques- 
tion is large or small because the corresponding rate constant is large 
or small. If w; is large because [C;] is large, 1) is true only if C; 
binds no catalytic component; otherwise the theorem of I applies; 
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viz., 3) the expression v, reduces to that for the steps 7 + 1 through 
7 +r, where the catalyst is returned to the free state in the 
(j + r)th step, and in this expression the concentration of bound 
catalyst is set equal to that of total catalyst. 4) If w., = 0 so that the 
overall process is irreversible, and w; is small, it can be shown di- 
rectly from (3) of I that 


j-1 
v, =i pe ) Gin (81) 


where G; is as defined in I. 5) If w_; is very small and w, is also 
small, so that the jth step is slow and irreversible, then from (66) 
or from (8) of I, 


it Wi; 1 
ysl ) ;——,) (82) 
Vv ( W; Ay 


k=-1 Wx 
which is valid whether or not w_, = 0. 


It should be obvious that, mathematically and hence physically, 
all of the above cases are not quite independent. The above classifi- 
cation is best for purposes of clarity and convenience. 

Setting aside, for the moment, the rather special case 5), it may 
be stated that in general in the above cases the rate expression has 
not reduced to that of a single limiting step. For example, in the case 
of the slow jth step, case 4) expressed by (81), the overall rate is 
indeed that of this step, it is in fact equal to the rate of any other 
step so long as steady state conditions are maintained; however, the 
rate expression has not reduced to that of a single step. According- 
ly, in the above situations, In v, vs. 1/T need not give a simple Arrhe- 
nius plot and if it appears to do so (and the following statement in- 
cludes 5)) the T-dependence is still not governed by a single A H# 
characteristic of a particular step. The rate expression does reduce 
to that for a single step in these cases, derivable from those above: 
a) If w; is large and 7 = 2 the first 7 — 1 steps amount to just the 
first step. b) If w_; is small and 7 = 1 the rate expression for the 
first 7 steps is again that for the first step alone. c) If [C;] is large 
and C; binds a catalyst which is regenerated in the (7 + 1)st step, 
i.e., 7 = 0 in statement 3) above, the rate expression is that for the 
(j ‘+ 1)st step. d) If the slow, irreversible step is the first one this 
is a degenerate case of (82) which gives v, = w,, the same result as 
b) above. The situation may now be summed up: In the general case 
the T-dependence of v, is given by (77). The characteristics of the 
In v, vs. 1/T curve is a problem for further work for the denomina- 
tors of both terms in (77) are strongly T-dependent. For the case 
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4), expressed by (81), 6 In v,/d 1/T will be similar to (77) since the 
derivative of the first factor in (81) will produce A H°,;. and Gin 
is of the same form as (66). For the case 5) it is evident by inspec- 
tion from (82) and (76) that 


Olnv, 


=== AH, i= jt ° 
Sp (AM + AHY)/R (83) 


Cases a) through d) require little discussion: When the conditions 
of these cases are met v, reduces to the rate expression for a single 
step and the T-dependence of v, is governed by the AH? of this 
particular step, provided the concentration terms are temperature in- 
dependent (see discussion to follow). 

The physical interpretation of (81), (82) and (83) is simple to 
discover: G;,, in (81) is, according to (7) of I, the function which 
multiplied by [X*;.] gives v,; the factor multiplying Gj» in (81) 
is, by (5) of I, the equilibrium concentration of X;.. Equation (81) 
then gives v; as the rate of steps 7 through n when Xj. is supplied 
to jth step at its equilibrium concentration. Equation (82) similarly 
gives v, as the rate of steps 7 through n when X;., is at its equilibrium 
concentration, but under the conditions in 5), this reduces to the rate 
expression for the jth step, by 2) above, and w; appears instead of 
a general term such as G;,,. Actually then, (82) is the rate expres- 
sion for the slow, irreversible jth step, but for reasons immediately 
to become obvious has been kept separate from a)—d). 

It is a point of some importance that in the situation 5), ex- 
pressed by (82) and (83), the plot of In v, vs. 1/T will be, to the 
_ extent that 4 H°,;.. and A H;? are temperature independent, a straight 
line. In some quarters, as already mentioned, the slope of this line 
would be interpreted as a measure of 4 H,;*. The error so committed 
is clear from (83) ; the slope is determined not by A H;* alone but by 
A H,? plus the sum of all the standard heats of reaction, A H;,° , from 
1 through 7 — 1. The physical situation is obvious: The condition 
for a “limiting step’? has been met, in fact the rate of the slow, irre- 
versible step is such that all preceding steps are essentially at equi- 
librium and the overall rate, (82), as already discussed, is just that 
of the slow step but this is not sufficient for the T-dependence of v, 
to be governed by A H,;? alone, for the equilibrium concentration of 
X;. is plainly temperature dependent. This is but an aspect of a 
more general problem already mentioned. When the rate of an over- 
all process is given by the expression for the rate of a single-step 
irreversible process the expression is the product of a rate constant 
and one or more concentration terms. These concentration terms 
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may themselves be temperature dependent: This may result from the 
temperature dependent distribution of any reactant or catalyst among 
active and inactive forms or among forms of varying activity; this 
includes inhibitions of all kinds. From results given in I, discussion 
p. 47, it follows that such distributions may be regarded as equilibria, 
though the main sequence of reactions are in the steady state. Al- 
though this had not previously been demonstrated, the assumption 
of equilibrium for such distributions has been the usual one in the 
past (Morales, 1947 and references therein, references in I). The 
T-dependence of such concentration terms thus enters through the 
equilibrium constants which describe the distributions mentioned 
above. More important, and this feature threatens the validity of 
many interpretations of “limiting steps and temperature coefficients” 
in in vivo systems, in open systems, every concentration is tempera- 
ture dependent. Each concentration in the system is maintained by 
outside or environmental “sources” and “sinks.” The temperature 
dependence in these cases is complex and enters through the diffu- 
sion factors of the environment and the system, permeability factors 
and the rate of the overall process in question (cf. examples in II 
and the many examples in Rashevsky, 1948). 

The situation treated by M. F. Morales (1947) is a special case 
of c) above, the speciality resulting from the fact that in his case the 
compound which binds the catalyst is X;. and in order that [X;-.] 
be large it is necessary, in our terms, that R, be small. His case may 
also be viewed as a special case of (82) wherein the expression for 
the equilibrium value of [X;_,] is obtained by solving the appropriate 
mass-action problem subject to the constraint of the material bal- 
ance on the catalyst, equation (4) of I, and the condition that R,, be 
small. Morales carried the problem further than we have here indi- 
cated by considering the effective total concentration to be T-depend- 
ent by virtue of reversible denaturation. It should be noted that if 
the step in which the catalyst is regenerated is neither slow nor irre- 
versible, statement c) still holds; in fact two such instances, equation 
(22) of I and discussion p. 39 of I, have already been given. 

Once again we invite attention to the fact that by and large the 
“limiting conditions” discussed above can be met, in general, by con- 
centration arrangement, without restricting the magnitude of rate 
constants which are inherent properties of the system. With this 
goes the important, if unfortunate, fact that in dealing with in vivo 
systems such arrangements are rarely possible. In the strict biologi- 
cal context the concentrations involved are, for a normal system un- 
der normal environmental conditions, be they internal and/or exter- 
nal, as much an innate property of the system as are the rate con- 
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stants. Finally, if one could always deal with the T-dependence of 
rate constants most of the problems discussed above vanish. In gen- 
eral however the biologist must deal with rates, and these under the 
conditions existing in the system. 


Some General Remarks. 


There are many qualitative concepts which can be discussed in 
terms of the formulation given here. A number of these have been 
introduced by A. C. Burton (1939). For example Burton speaks of 
the “cost of transition’ between steady states. In terms of “heat” 
this quantity has already been indicated in (60), or, in terms of free 
energy, and this seems more appropriate to “cost,” in (47). Similarly 
the quantity F, may be regarded as the “cost of maintenance,” it is 
the free energy expenditure per unit time associated with the steady 
state and will be indicated in IV following, to be the rate of dissipa- 
tion in the system. Burton further points out that the heat capacity 
of a steady state system may be expected to involve the “‘heats” asso- 
ciated with the temperature-induced concentration changes. This is 
of course true but is not a unique property of the steady state for 
such results are well known for equilibrium systems (Zemansky, 1943 
p. 351). The computation is easily indicated for the system discussed 
in this paper, for 


dH=> md He +S. Hid m (84) 
k k 


where the mole numbers 7; have been defined and the H; are the par- 
tial molal enthalpies: If the variation (84) is due to a temperature 
change dT and, see discussion following (20), the changes in the 
mole numbers of the C; and B, are neglected and only the [X;*] con- 
sidered to change, (84) yields 
aH\ _ a ~ / [Xe] 

(— ik Cs Sm Co + VE H,( 5 T F 5 (85) 
‘The first summation, in which Cp, are the partial molal heat capac- 
ities, is the heat capacity of the system as if it were an inert mixture, 
the last summation is the contribution due to the change in the X-con- 
centration: The 0[X;°]/eT are readily obtained from (5) of I writ- 
ten as 


Kix 


[Xe] = 


k 


ry] Vis (86) 
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where V.;,, is defined in I, and using (68) and (69) of this paper. It 
is not worthwhile to carry this out here. 

There has appeared in this paper no quantity which corresponds 
to what Burton calls the “potential energy.” The “potential energy” 
for an equilibrium system is the thermodynamic potential or F’, which 
is a minimum at equilibrium: If such a system be displaced from 
equilibrium and abandoned, F will decrease until the minimum value 
again obtains. If a steady state system be displaced it will also re- 
turn, so far as the [X;°] are concerned to the original condition; the 
system discussed here does not return to its original F-value; the 
difference between the F-flux at any stage of the recovery and that 
in the steady state, (26), has been discussed but the difference in the 
F-values has no real meaning in the sense that it does for an equilib- 
rium system. It is a fundamental problem to inquire whether there 
exists for a steady state system a function which has the property 
(which the various thermodynamic functions have for an equilibrium 
system) that it is a maximum or minimum in the steady state. This 
question, to the author’s knowledge, first raised by A. J. Lotka (1925 
p. 821) has been discussed in interesting and definite terms by F. 
Zwicky (1933). The specification of a steady state system is a dif- 
ferent problem from the specification of an equilibrium system in that, 
for an n-step system, 2n w’s must be prescribed; this involves 2n 
rate constants and as many concentration terms as are included in 
the w’s. The condition for the steady state is, (23), vx = vi for all k 
and there are elaborate tests (Lotka, 1925; Garder and Barnes, 1942 
and refs. therein) for the stability of the steady state but it remains 
a problem to relate these conditions to any energetic criteria. 

We have discussed only “straight-chain” sequences of reactions 
in this paper. The extension of these formulations to the branched 
or coupled systems treated in I offers no difficulty, for example F(t) 
is still given by summing A F; v; over all steps, cf. (8), (16). In 
connection with this type system and closely related to the conclusion 
from (21) and (22) is the following problem: Consider a system in 
which a substrate is being oxidized with the concomitant formation 
of an “energy-rich” compound and which is at least triply branched 
and coupled in that 1) the reduced coenzyme is reoxidized via a hy- 
drogen transport chain, cf. I; the energy rich compound donates a 
group, usually a high-energy phosphate, to an acceptor; 2) the or- 
ganic moiety, divested of the donated group, suffers some further 
fate via a chain of reactions and 3) the now energy rich acceptor, 
X , can possibly divest itself of the accepted group via one or more 
pathways which may include wasteful splitting or spontaneous de- 
composition. Now the steady state (or time dependent) concentra- 
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tion of X is determined by the rates of all of these processes. A given 
process which requires or must be “driven by” X will not take place 
unless the concentration of X is sufficiently great; in fact this process 
will, according to (22), run backwards. Thus if the rate of utiliza- 
tion of the substrate is not sufficiently rapid the oxidative process 
cannot “drive” the process in question; specifically the rates of the 
primary dehydrogenation and process 2) above must be rapid, mean- 
ing that, in the steady state, the forward w’s of these steps must be 
great. The situation is simply that the concentration of X is deter- 
mined by the rates involved and that free energy, residing in the sub- 
strate and ultimately available in the form of energy-rich acceptor 
may be sufficient to “drive” the process in question; this does not oc- 
cur if the chemical potential, given in its dependency on concentration 
by (2), does not become sufficiently great. The acceptor is not nec- 
essary to the above picture and frequently the energy-rich compound, 
such as an acyl-phosphate, per se participates in synthesis etc.; the 
general argument remains the same. On the above basis a given syn- 
thesis, dependent upon an oxidative process, will not occur unless rate 
of the oxidative process is sufficient, in fact by (22), degradation will 
occur. Further there will be competition, for “energy,” in that rapid 
syntheses keep the X-concentration at a level insufficient for some 
slower processes to occur. Likewise there will be competition between 
“maintenance” and function, to the extent (which is not known) that 
function depends upon the concentrations of energy-rich intermedi- 
ates in the same way as do rate processes; there may also be compe- 
tition between any metabolically supported processes subserved by 
the same group of energy rich intermediates. It is a simple problem 
to verify the above statements in principle. Their experimental im- 
portance is, of course, another question; a complete discussion can- 
not be given here. A presumptive case of the failure of an oxidative 
process to support function due to insufficient rate is known (Tschirgi 
et al, 1949). An example of competition between syntheses is that 
between the syntheses of virus and normal protein. It would be in- 
teresting to investigate the energetics of development (Tyler, 1942) 
and similar processes from the above point of view. 

It will be mentioned here that the formulations in this paper are 
clearly applications of classical thermodynamics to systems in which 
changes, irreversible and of any rapidity, are occurring. It may seem 
unnecessary to make such a remark at all, but the idea is sufficiently 
prevalent that the methods of thermodynamics are not available for 
non-equilibrium systems and irreversible processes and the failure 
to recognize thermodynamic methods in unfamiliar context suffici- 
ently frequent that this will be discussed at the close of IV following. 
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The investigation initiated in III is continued. The irreversible pro- 
duction of entropy and the rates of dissipation in diffusion fields and due 
to chemical reactions are discussed. The rates of change of free energy, 
entropy, etc. are indicated for type-systems of biological interest in 
which composition and reaction rates are a function of position. Some 
consequences of Onsager’s generalization of Fick’s Law are discussed in 
terms of maintenance of stationary, non-equilibrium concentration dis- 
tributions, transport of a solute against a concentration gradient, and the 
dependence of these phenomena upon metabolism. 


In the introduction of part III (Hearon, 1950) it was indicated 
that the dissipative or irreversible aspects of processes are of some 
general interest for the reasons there stated and frequent reference 
was made to “dissipation” in steady state systems. It was further 
proposed there to discuss, in this part, open systems and systems in 
which space gradients of the thermodynamic quantities occur. We 
begin here by discussing the irreversible production of entropy. 


1. The Irreversible Entropy and Dissipation. 


In principle there will be nothing new in this section: There are, 
as has been indicated in III, many excellent treatments of what is es- 
sentially this same topic. Th. De Donder first introduced explicitly 
the “uncompensated heat” of a process (1936 and references therein). 
He did not employ the term irreversible entropy but related the un- 
compensated heat to what he termed the affinity (see also Van Ryssel- 
berghe, 1935). The affinity is shown to be 


d * 
Affinity = ds 
oF ") 
=-(> ) =— Bom, 


where q’ is the uncompensated heat, & is the extent of the reaction and 
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the members of the second equality have been defined and discussed 
in III [equation (5)]. The development leading to the result (1), 
while an essential part of the logical apparatus in the remarkable De 
Donder-Van Rysselberghe exposition of thermodynamics, is not very 
direct. Furthermore the formulations of these authors are restricted 
to systems the composition of which can be specified by one or more 
variables &;; in particular they are not readily adapted to processes 
such as diffusion. In a general treatment of irreversible processes L. 
Onsager (1931a, b, 1932) constructs for transport processes, notably 
diffusion (see also Harned, 1947), a dissipation-function: Twice the 
dissipation-function gives the time rate of change of the irreversible 
entropy. In Onsager’s treatment there is no necessity for the dis- 
tinction, which will be made in what follows, between total and irre- 
versible entropy. 

In the general treatment of irreversible processes by C. Eckart 
(1940a, b) the separation of the irreversible entropy change from the 
total entropy change is actually made but not quite in these terms. 
It is perhaps permissible to state here that Eckart’s very general ob- 
ject (1940b) was to inquire for the restrictions imposed by the sec- 
ond law of thermodynamics upon the equations, which he calls the 
thermodynamic equations of motion (cf. Onsager, 1931la, b), which de- 
scribe irreversible processes} and to derive equations of the form 
[Eckart’s equation (1) in the notation of this paper], 


a i sav+ J at) -de= J eav, (2) 


which, with the constraints of the second law of thermodynamics, 
give 
G20, (3) 


where S is _the total entropy per unit volume, g the vector of heat 


flow and d o the vector element of surface with direction of the out- 
ward normal. The quantity G in (2) and (3), given explicitly in 
Eckart’s equation (26), may be identified at once with what we, fol- 
lowing R. C. Tolman and P. C. Fine (1948), will call the (time) rate 
of irreversible production of entropy per unit volume. J. Meixner 
(1941, 1942, 1943) and B. Leaf (1946) have presented treatments 
similar to that of Eckart (1940b). Both of these workers use the term 
dissipation in the sense that we will immediately employ it. 

In view of all of this it merits some comment that what we hope 


+We have already referred, in III, to his conclusion that these restrictions. 
are fulfilled by the formal kinetic rate equations for chemical reaction. 
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to achieve in this section and in the remainder of this paper is some 
sort of synthesis of the above-mentioned results in a fairly direct, if 
simple, treatment. This will be carried out in a manner strongly di- 
rected toward our interests here, which are not too general, and 
with regard for rendering complete certain statements and formu- 
lations in III. More especially there will be discussed and employed a 
simple expression for the irreversible production of entropy which 
applies to a large class of processes. It includes, for example, many 
of the individual cases treated by Tolman and Fine (1948). It should 
be made clear that generality was no primary concern in their work 
the burden of which was to compute from the familiar integral 
AS =f{dq/T, taken over a reversible path, the irreversible production 
of entropy for a host of individual cases, to provide a logically satis- 
factory basis for this procedure and to examine the rigor with which 
the rate of irreversible production of entropy can be computed. 

Since most of the above cited investigators have made little ref- 
erence to previous and related work, the relations which are briefly 
and incidentally pointed out in what follows may be helpful outside 
of the immediate context of III and of this paper. 

Although the total entropy, S, is an extensive quantity and is 
in this sense an additive property (see Epstein, 1937 for a careful 
statement) it is, according to the second law of thermodynamics, not 
a quantity which is conserved in irreversible processes. Thus in con- 
sidering the total entropy change of a system it is necessary to ac- 
count for the entropy reversibly communicated to the system across 
the boundary and the entropy produced within the system. Denote 
by S, the reversibly transported entropy and by S* the entropy gen- 
erated within the system by irreversible processes. The total entropy 
change is then (cf. La Mer et al, 1949) 


dS=dS,+ds". (4) 


Let dq be the thermal energy received by the system in any varia- 
tion. Let dq* be the heat irreversibly evolved within the system: this 
is the uncompensated heat already mentioned. The second law states 
that 


TdS = dq. (5) 
The equality statement of the second law (De Donder, 1936) is that 
T dS =dq + dq’. (6) 

From (5) and (6) we get 
dq*=0. (7) 


For a reversible process the equality in (5) and (7) holds. If no mat- 
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ter is transported into the system and entropy is carried into the sys- 
tem by heat transfers alone, then we find 


dS, =dq/T ; 
i= de. 


The last-equation in (8) obtains whether or not matter is exchanged 
with the environment. If the quantities in the last equation of (8) 
be referred to a unit volume of system and the time derivatives formed 
from the differentials there is obtained Eckart’s quantity G [equa- 
tion (1)]. The inequality (7), demanded by the conventional second 
law (5), is equivalent to (2). 

It should be mentioned here that we have (1949a, b, referred to 
as I and II), in purely kinetic discussions and in III, freely used the 
term irreversible to describe a reaction in which the reverse or back- 
ward rates can be neglected: This is not to be confused with thermo- 
dynamic irreversibility, here discussed, to which it bears little rela- 
tion for all reactions under. non-equilibrium conditions are thermo- 
dynamically irreversible. We will also assume as true the coroilary 
of the second law that processes which are proceeding at a finite rate 
are irreversible but note in passing the hypothetical exceptions dis- 
cussed by J. Meixner (1942) and the relativistic situations discussed 
by R. C. Tolman (1934, p.134, 319). These exceptions need not con- 
cern us here. 

Now the method of computing dS* or, for a finite change, 4S* 
will naturally be determined to a large extent by the system and pro- 
cesses in question. In general AS* will consist of terms giving the 
contribution from diffusion, spontaneous chemical reaction, mix- 
ing, heat flow, work against viscous forces, etc. As has been shown 
(Tolman and Fine, 1948) the calculation of these terms can usually 
be made by computing /{dq/T over a reversible path connecting the 
initial and final states of the system. One has only to invent a process 
entailing the reversible absorption of heat which will produce the 
same change in the system as that produced by the irreversible pro- 
cess in question: The limitations imposed on the use of this method 
have been carefully discussed (Tolman and Fine, 1948). In the case 
of thermal effects alone this procedure is simple. For other processes 
greater resource is required on the part of the computer to conceive 
of a process involving reversible heat exchange which will produce 
the desired effect. The following method is most frequently avail- 
able; in particular it is of use for the processes of diffusion and 
chemical reaction in which we are chiefly interested. 


Consider a system, in which any process is occurring, to be 


(8) 
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placed in thermal contact with a huge, infinite, if necessary, heat res- 
ervoir. The temperature, 7, of the system and reservoir will be the 
Same and the combination, system plus reservoir, will operate at con- 
stant energy E and constant volume V. It is clear that if the change 
in entropy be computed for the combination, taking account of the 
fact that the reservoir can lose energy and volume to the system, the 
entropy change so computed will be just the irreversible production 
of entropy in the system. Consider the entropy of the system, 
S.(T, P, xi), to be a function of the temperature, pressure and 
any number of parameters, x; , which specify the composition. Con- 
sider the entropy, S,(£,, Vz), of the reservoir to be a function of 
the volume and energy. We exclude for the moment exchange of mat- 
ter between the system and reservoir. A change in the parameters 
x; produces a change, 


0 Sa 


i ON 


ax: , (9) 


in the entropy of the system and changes, 


OK, 
Cs, dx;, 
ane 
(10) 
oVa 
dV.e=> ax;, 
pe ae 


in the energy and volume. Because of the constancy of total energy, 
E=E,+ E,, and volume, V= V, + Vz, the changes (10) cause 
equal and opposite changes in the reservoir, 


dE >= dH “5 
(11) 
dV, =>=—_ ave 
with the result that the entropy of the reservoir changes by the 
amount, 


as aici 
ds dB, ti OV 
0 Ks 0 Vs, 
(12) 
1 
= diate dV a: 


For the total entropy change, dS, + dS,, there is obtained, by using 
(10) and (11) in (12) and adding the result to (9), 
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0 Sa p bet oe PONV a 
,= > — dx, - — dz,—— dx;. (13) 
But at constant 7 and P, 
aF, 
asa (dE, + PdVa)/T——— (14) 
and (13) gives 
dF.) 2, 
(dS) py =— SP = as (15) 


The above is in all essentials the same as the method employed by J. 
Mayer and M. G. Mayer (1940) to show from the second law state- 
ment, (dS)z, = 0, that F must decrease in spontaneous processes. 
In the above derivation it was assumed that 7 and P were uniform 
throughout the system. The composition was not assumed uniform 
and the irreversible processes taking place in the system might in- 
clude mixing, redistribution of concentration gradients with trans- 
port from regions of higher to regions of lower chemical potential, 
etc.t 

Consider now the situation that the reservoir contains all of the 
chemical species contained in the system and that exchange of mat- 
ter may take place between the system and reservoir. Let S,, EH, and 
V, be functions of 7, P and the mole numbers ;, and S; be a function 
of E,, V, and n;,. In place of (9) we now have 


dSc=D Sind Nia, (16) 


where the Sin are the partial molal entropies in the system and (12) 
is replaced by 


dS) = By + FAV) — FS tn Nw, (17) 
where the w;, are the chemical potentials in the reservoir: 
0 Sp 
( al Ie =— uin/T. (18) 


From (17), (11) and (16) there results 


tIt may then be true that (9) and (10) may not be expressible by a finite 
number of «; and a space integral may be necessary but only the changes dS,, 
dE, and dV,, whatever their form, appear in (18) and the argument holds. 
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Lid 


it P 
(dS) , =dS,—— dH,——d = i 

E,V T T Vs ~ T dni» 

(dF, + dF,) np ee 
= : 

Since in (19), dS, is given by (16) and dE, and dV, by 

dE. =D Bia Unin (20) 
and 

dVa = > Ve Nias (21) 


where EH, and V;, are the partial molal energies and volumes, and the 
relations, 


tia = Hie —T Sia, 
is za ¢ (22) 
A,=Enwt PVia 


obtain among the partial molal quantities in question, it is readily 
seen that the contribution to dF, + dF; due to transport of matter 
from the reservoir to the system will vanish if wic = min for alli. Thus 
if the chemical potentials are equal in the system and reservoir equa- 
tion (19) still gives the change due to spontaneous processes in the 
system. If the composition of the system is not uniform, (16) must be 
computed by considering the dn;, to be the change in the mole numbers 
for a unit volume and integrating over the volume of the system. The 
point of importance is that if the surface values of ui, and ui, are 
equal, (19) still gives the change due to irreversible processes in the 
system since that due to transport across the surface will vanish as 
discussed above and (19) gives the change due to transport, in the 
system, along gradients of chemical potential, chemical reaction, etc. 
This case of the reservoir containing matter which can be exchanged 
with the system has been discussed because of the strict formal ana- 
logy, which later becomes clear, to the cell imbedded in an infinite 


medium: 
From (15) and from the familiar 


AF=AH—-TAS, (23) 


which holds at constant T, it is seen that AS and A S* are identical 
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if A H = 0, i.e. for a process for which the enthalpy change is zero 
the total entropy change and the irreversible production of entropy 
are the same. This is in fact the case for the mixing of isotopes or 
of perfect gases, for in these processes 


AS=—AF/T (24) 


(Mayer and Mayer, 1940). This will be shown in the next section to 
be the case also for diffusion in perfect, dilute solutions. Using (4) 
in (28) gives 


AF=AH—T(AS,+ AS"), (25) 
and comparison to (15) suggests that 
A A Ss (26) 


It may be argued that at constant T, when (23) and (25) apply, the 
heat A H is absorbed by the system from the environment reversibly. 
That this fact-is included in (15) may be easily shown in a manner 
which will exhibit the T-dependence of A S*. Differentiation of both 
sides of 


—TAS*=AF (27) 
with respect to T and using the relation 
OAF 
Fm esa me vai AiO") (28) 
show that 
CAR eure (29) 
iT sa eae 
But also 


O Ais ci eo Ad yal, 
oT a i 
Comparison of (29) and (30) gives (26). 

The quantity TdS*, or T A S*, gives, cf. (8), the uncompensated 
heat and represents the amount of dissipation associated with the 
process in that this amount of energy, formerly available for useful 
work, has become completely unavailable for work at the temperature 
T in question [cf. discussion by Zemansky, 1943, p. 167 wherein the 
amount of energy which becomes unavailabe is 7) A S(system) + 
T, AS (environment), where 7, is the temperature of the coldest avail- 
able heat sink]. This is easily seen (LaMer et al, 1949) by substitut- 
ing dq from (6) into the first law of thermodynamics (see also Tol- 


d. * 
man and Fine, 1949). Thus T ( = 


=AH/T?. (30) 


) gives the rate of dissipation: 
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this is a form of Kelvin’s principle of degradation of energy. 


2. Free Energy Change and Dissipation in Diffusion. 


Consider the diffusion of a solute from a point where the chemi- 
cal potential is 4(1) to a point where the chemical potential is (2). 
For one mole the change in free energy is u(2) — mw(1) and using 
(15) in the form of (27) 


PEAS ia) e ge (1 i) (31) 
Taking yu to be of the form [(2) of III] 
p= wtRTINC (32) 
gives 
Zz, 
Pl eT ye. (33) 
C(1) 


Since yw is related to the partial molal enthalpy H by 
Bp de o0 ue JL = 
—_— = =— H/T? 34 
oT Oel: / Pd 
and uw? is independent of concentration (cf. III), it is implicit in (32) 
that the transfer of solute from one concentration to another is not 
accompanied by the absorption of heat. Thus A S* in (88) is iden- 
tical with AS. 
Consider an infinitesimal cube located at the general position x, 
y, 2. The rate of transport across the face in the y-z plane in the 
positive x-direction is J, dydz, (Moles/time) , where J, is the x-com- 


ponent of the diffusion current J. The change in free energy, for one 


Ou 
mole, in being transported the distance dx is ae dz. Thus the time 
x 


rate of change of free energy per unit volume, when there is no ac- 
cumulation in the volume dx dy dz, due to flow in the positive x-di- 


0 
rection, is J, —-. Adding the similar terms for flow in the y and z 
z 


directions gives for the time rate of change of free energy per unit 
volume 

aF es 

aiinens geo od LO J, (35) 


where F is the free energy per unit volume. Thus the total rate of 
change due to diffusion is 
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d. F 2 > 
ee Vuod dV, (36) 
dt V 


where F is the total free energy. From (15) and (386) there is ob- 
tained 


d S* 

dt 
as the total rate of dissipation due to diffusion. It can easily be shown 
that (37) gives the work per unit time against frictional forces. Let 


=— [ vp-Jav (37) 
A: 


J be given by the simple Fick’s Law 


eae Saeccle (38) 


The velocity of diffusion is 
es DN CG V (39) 
oe a Sa eS he 
v L 


Now one may, following E. A. Guggenheim (1929), regard V u 
as the force on a mole of solute (cf. Hartley, 1931; Harned, 1947) 
and identify RT/D as the “frictional coefficient.” The work per unit 
time, per mole against frictional resistance is then 


TF i ee ees (40) 


Since there are C moles per unit volume the work per unit time per 
unit volume is 


yee nt (41) 


Integration of (41) over the volume in question gives (37). N. Ra- 
shevsky (1948) has derived (41) in the form DRT(VC)2/C with- 
out use of the chemical potential but by invoking the Stoke’s Law ex- 
pression for the force of resistance upon a spherical particle. The 


result (37) holds whatever the form of J so long as this is the cur- 


rent of diffusion. Thus J need not be given by (38), which was used 
only to show that (37) also represented degradation of energy by 
friction, but may be taken from Onsager’s (1932, 1945) generali- 
zation of Fick’s Law which we will later mention and use. The cur- 


rent J however must not include a convection term for then the 
change in free energy would include a term due to transport into the 
region which does not represent dissipation. When convection does 
occur (37) may be derived by considering an infinitesimal cube which 
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moves with the fluid and the volume V of integration in (37) would 
likewise be considered to move with the fluid. 


3. Free Energy Change and Dissipation for Chemical Reaction. 


With the aid of (15) and the formulations of III this topic may 
now be disposed of rather simply. From (15) 


GSH sn lide 
dt T dt 


(42) 


Lf Geen 
and i is in general to be taken from equation (8) or (14) of III 
or, if the steady state obtains, from equation (16), (20) or (43) of 
III. The total rate of dissipation, T ( 


* 


di ), due to chemical reac- 


pre Les dF : 
tion is — i from equations (8), (4), (16), (20) or (48) of III. 


Similarly equation (26) of III gives the difference between total rates 
of dissipation in the time dependent and steady states. Equations 
(44), (45) and (46) of III give the difference between the rates of 
dissipation in two steady states when the transition between steady 
states has been induced by one of the various means discussed in III. 
Equations (51) and (52) of III, with a change in sign, give the 
T-dependence of the rate of irreversible production of entropy and 


d. * 
the rate of dissipation respectively. The rate of dissipation T ( ) 


is required by the second law, equations (7) and (8), to be positive. 
By writing (42) in either of the forms in the notation of IIT 


T —F; 
dt 
(43) 
don 
r( 2S") = Fy, 
tei. 


where the second form holds for the steady state, and recalling that 
in III it has been shown that F < 0, F, < 0, one obtains from a 
slightly different approach confirmation of the fact, cf. discussion 
p. 63-65 of III, that the rate laws of formal kinetics are in accord 
with general thermodynamic principles (cf. Eckart, 1940b). 

It should be mentioned that LaMer et al (1949) have derived 
the equivalent of (42) for chemical reaction by employing their pow- 
erful and elegant rule of potentials (p. 618 LaMer et al, 1949) which 
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cannot be discussed here but to which we wish to direct attention. 
We return for the moment to the discussion of section 1 of this 
paper. Th. DeDonder (1936) gives dq*/dt as 
dq* = dq* dé 


: (44) 
dt dé dt 


Ki: dF 
From (1), cf. (5) and (6) of III, this is — <r DeDonder defines 


dq*/dt as the “power of the system” and shows that it must be posi- 
tive: It is readily seen from the discussion in section 1 of this paper 
and in III that this is equivalent to our demonstration that F < 0, 
FEO 

It is worthy of mention that the relation discussed above between 
dissipation and free energy change gives some insight into the much 
discussed (Lipman, 1941; Barron, 1943) situation that biological sys- 
tems are usually characterized by processes which take place in many 
‘nearly reversible” steps and which are thus less wasteful. Simi- 
larly the conservation of free energy resulting from the coupling of 
processes such as oxidation and phosphorylation (Lipman, 1941; cf. 
III) receives a somewhat clearer interpretation in terms of dissipa- 
tion. 

In order to extend the results of this section and those of III to 
systems in which the concentrations of the various components are 
functions of position we give the following discussion: Consider an 
element of volume A V sufficiently small that the concentration may 
be considered uniform within the element. Let F be the free energy 
and 7, the number of moles of the kth solute in the element. Then 


F's Me; (45) 
k 


where the , are uniform in the element. If the volume element ap- 
proaches zero 


_— 


F i 
1 ee SS : 
Mer AV ~ bx Cx, (46) 
where F has already been defined as the free energy per unit volume 
and the concentration C; is the number of moles per unit volume, at 
the point in question, of the kth solute. In a system with space dis- 
tributions of concentration it is no longer true that the total free 
energy F' is a homogeneous function of the mole numbers of degree 
one; for the number of moles of the kth solute is 
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hy = | C.dV 
V 
and the total free energy is 
LO | aR S me C, dV. (47) 
V VK 
Similar statements apply to the other so-called extensive properties. 
The free energy per unit volume from (46) is clearly a homogeneous 


function, of degree one, of the C;, and the usual theorems based on 
this fact apply. In particular we have 


0 Bj 0 bk 
SKE: = == @;, 48 
: "OG : SaC) (48) 
and 
> C.d w= 035 dT =4dP—0; (49) 
k 


which are, for a unit volume of system, forms of the familiar Gibbs- 
Duhem rule (cf. Vegard, 1910). Thus from (46) and (49) the rate 
of change of free energy per unit volume, at constant T and P is 

7 d Cy 

fe Ser 
in which «4, and C; are evaluated at the point in question and dC;/dt 
is due solely to chemical reaction. If dC,,/dt be evaluated in terms 
of the kinetic rate equations the result is obviously equation (8) or 
(14) of III with the factor V deleted. In fact it was indicated in III 
that AF;v; gives the rate of change of free energy per unit volume 
due to the ith step. The factor V appears in (8) of III because of 
this. If (50) so expressed be integrated over the volume of the sys- 
tem the total rate of change of free energy due to chemical reaction 
is obtained. If the steady state obtains (50) is given by equation 
(16), (20) or (48) of III with the factor V deleted. 


(50) 


4, Metabolizing Systems. 


Consider a system in which any number of substances are being 
produced by chemical reactions at the rates, per unit volume, q; 
which may or may not be independent. The equation of continuity is 


re ac 
EDA era (51) 
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And the total rate of eo of free energy is 


ak F_ 0 Cy 
FaV= fe —dV, (52) 
a dt alt ere 

where 0C;/ot is due, by (51), to diffusion and reaction. Combining 
(51) and (52) and applying the divergence theorem gives 


<= | Smut [ SVmndV— [ Satedo, (63) 
Vk v & o &k 

where oa is the surface enclosing the volume V. Since q < 0 for con- 
sumption the first integral is a contribution from chemical reaction; 
the intergrand is given by (50) and may be expressed in terms of 
the formulae of III as discussed in section 38. The second integral is 
by (36) the contribution due to diffusion. The surface integral is 


the contribution due to influx of matter at the rate ee do moles per 
unit time through the element of area dw, each mole of the kth so- 
lute bringing in free energy in the amount u,. It is clear that a simi- 
lar procedure will give the rate of change of entropy when 


dt “at Se 
is used and so on for any thermodynamic quantity, for the argument 
leading to (46) may be applied to any quantity and the corresponding 
partial molal quantities. 

What can be done to carry the problem further depends entirely 
upon the solutions of the set (51). As is well known, if the gq, are re- 
lated by the fact that any q, depends upon other concentrations C; , 
even the steady state solutions of (51) are usually unavailable. When 
Gd, is just proportional to C;, the steady state solutions are rather un- 


wieldy (Rashevsky, 1948). When the J; are not independent, as will 
be true when Onsager’s generalized law (1945) is used, situation 
arise, to be mentioned later, not previously treated. There are a few 
points however which are worth discussing with the aid of very sim- 
ple examples. 

It is clear from (52) that when the steady state has been attained, 
0 Cy 


dF 
er —=() fOr all K, =H is zero. The integrals in (53) extended over 


the volume of the cell, (system), and its boundary surface o must can- 
cel: the contributions from chemical reaction, diffusion along gradi- 
ents of chemical potential and influx of matter must balance. When 
one considers the medium, for which the equation of continuity is 


BS ease 
V k 
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m5 0:Cx 
=F (55) 
in place of (53) there is 
dF = =a 
(=) = [ 37 m-hav— Sere de (56) 
dt medium vieo&k o &k 


Equation (56) must also vanish when the steady state obtains. Here 
we encounter a difficulty pointed up in III: If among the boundary 
conditions on (55) is the usual condition that C; > C;° at infinity then 
it is considered that material is added to or withdrawn from the me- 
dium, by reactions in the cell, without changing the concentration dis- 
tribution. Thus there results the situation that material flows into 
the cell, reacts, products flow out; all with no change in free energy. 
In the system discussed in III the free energy change accompanying 
the conversion of reactants to products was admitted but, in the rate 
equations, changes in concentration were not admitted. In the prob- 
lem discussed here the situation is not so simple to remedy. The prob- 
lems associated with the “infinite medium” were early discussed by 
N. Rashevsky (1948, preface to first edition). 

Consider the simple case of a spherical cell of radius 7) in which 
the reaction C = B occurs. Let there be no gradients except in the 
membrane. The rate of change of free energy for the cell and envi- 
ronment due to diffusion of a single solute is, 


To+d 
aC. dln C 
(;) Sor (22) ‘dr, (57) 
dt /aitt r Or /m or 


0 


where (0C;/07) m is the gradient in the membrane which is assumed 
constant, D;., the diffusion coefficient in the membrane and 6 the 
thickness of the membrane. Equation (57) gives approximately 


Cz (7% +6 
(3) ee ro( | ree ey poten (58) 
dt / ate Or /m Ci. (10) 


Since C(7)) and C(” + 6) correspond to C, and C;, the external 
and internal concentrations, (58) gives 


Co 
Ce 
where hA(C, — C;) isthe familiar expression (Rashevsky, 1948) for 


the flow into the cell with h the permeability constant and A the area 
of the cell. Let the reaction be proceeding at the steady state rate 


(=) =— I A(Co— Ci) RT In (59) 
dt diff 
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v, = k,[C] — k.[B] and let [C]. and [B], denote the external con- 
centrations. Now, cf. III, 


( a y Z v( ar + RT In a v, (60) 

and 

(F va = heA ((Cle— (€]) RT na 
[B]o 


h;A ((B] — [B]o) RT In ——. 
+ hgA ([B] — [B]o) 1B 


(61) 


In the steady state 


Vv, = hcA ([C]o— [C]) = heA ([B] — [B]0)-. (62) 
Thus the sum of (60) and (61) gives 


( = ) ces = ( AF® + RT In rao) Vs; (63) 


where AF” is defined in III. Equation (63) gives the rate of change 
of free energy for the conversion of C to B at the external concen- 
trations [C],. and [B], at the rate v,. The physical situation behind 
(63) is certainly clear. This simple example suggests at once that the 
method of avoiding the difficulties caused by the “infinite medium” 
(which is a necessity for solving the boundary value problems of 
diffusion in closed form) is to employ the approximation method of 
Rashevsky (1948). This will involve certain problems but, in gen- 
eral, one has only to set up approximate gradients and find approxi- 
mate volumes to replace the volume integrals. 

Thus far we have wherever possible left the form of the diffu- 
sion currents J; unspecified. Where the form has been given the 
conventional Fick’s Law was assumed. We wish to adapt the perhaps 


unsatisfactory procedure of drawing heavily upon results not yet pub- 
lished (Hearon, 1950b) in order to discuss a most interesting case. 


If one, following Onsager (1932, 1945), takes the J; to be of the form 
Jx => Diy V C; (64) 

j 
where D;; is the conventional diffusion coefficient and D;;, i #7, are 


the “drag coefficients,” the equation of continuity becomes, for the 
steady state, in the cell 
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> Di? V? Ci§ =— qe, k=1,2,----. (65) 
j 


It has been assumed that the D;;' are essentially independent of con- 
centration. Consider a spherical cell of radius 7 in which the q% 
are constant. The set (65) gives for the internal concentrations, C;' 


Qk 7? 
CS San + Ay, K=1,2,---, (66) 


where r is radius vector from the cell center, the A, constants and 
the Q; are for the case of two solutes, 


— Q1= (— @ Das’ + Qe Diz) / (Dur* Doz’ — Dye‘ Dz;') ; 
— Qs = (— @2 Du? + Gs Dor) /(Dur‘ Doo — Diz’ Dar‘). 
The external concentrations C;° are given by 
Ci? = C,° + Bi/r, (68) 


where C;° is C,° at r= o and 


(67) 


T° 


B= 

‘ (69) 
7 oe 
2 2 3 ’ 


where Q;° is the same function of the external coefficients D;;° that 
Q; is of the D;;'. The A; in (66) are given by 


To” 
3 


where Q,” is same function of the h,; that Q; is of the D,;* and the 
hy; are generalized permeabilities having the same significance for 
flow across the membrane that the D,;' have for flow in the cell. For 
example the flow from the cell across the membrane of the first sub- 
stance is 


To 1” 
Ay = Cy? + Qe” a + Qk ra HO; (70) 


Ia (C,' a Ci*) o te his (C.' = Cs") 0% (71) 


where the subscript zero denotes that the values are at r= 1%. As 
is obtained from A, by interchanging the subscripts 1 and 2 every- 
where in (70). Now let q, = 0 and q, be negative, corresponding to 
consumption in the cell of the second substance. Then 


Q,:=— Dt /(D3' Dy — Di" Das); 


: ; (72) 
Q.= Dis CDi" Da! = D323 Dz,'), 
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and it can be shown from (66), (68), (69), (70) and (72) that the 
concentration of the first solute, C,’, though this solute is neither con- 
sumed nor produced, is everywhere higher in the cell than outside. 
It is in fact readily verified from (66) with (68), (67) and (70) 
that C,‘(r) is a paraboloid with its apex at the center of the cell. 
The distribution is formally identical with that (Rashevsky, 1948) 
of a solute which is being produced in the cell. The distribution C,' 
depends upon Q, from (72) in such a manner that the difference 
between C,'(0) and C,‘ (7) depends upon gq, the rate of consumption 
of the second solute: If |q.| is increased the distribution changes so 
that C,‘ is even higher relative to C,° and the first solute flows into the 
cell against a concentration gradient. If the initial conditions are that 
att =0,C,' = C,° the “free energy of concentration” is given by, cf. 
(47), 


eae 4 
A F dene. =A4 af haat Cy’ T2207 — Seal! T° Mio c;° ) (73) 


oO 


where jf. is the chemical potential of the first solute at the concen- 
tration C,°; uw, and C,* are to be expressed as a function of r from 
(82) and (66). Equation (73) gives the change in free energy of 
the cell accompanying the achievement of the steady state distribu- 
tion from the initial distribution of the first solute. Similarly for two 
steady state distributions, 


A t ets =4an | Mie Cw r? dr— An | Maa Crs v2 dr (74) 


gives the “free energy of concentration” due to the first solute where 
Mx, and {42 are the chemical potentials at any point r in the cell for 
the first and second steady state, a similar notation is used for C,:‘ 
and C,,.'. Equations (73) and (74) give also the work which must 
be done upon the cell to effect such a change in distribution: It is 
easily seen that this is derived from free energy which would have 
been dissipated in the diffusion field of the second solute had the first 
solute not been present. From (66) 


V Ci r/r=—Q, 7/3; 

= (75) 
V C,'-r/r=—Q.7/3. 

Taking Q, and Q. from (72) and ip from (64) shows that 


Jj=J"=0. (76) 
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Thus there is no net flow of the first solute and hence no dissipation 
in its diffusion field, though VC, is everywhere different from zero 
except at r = oo. From (36) the rate of change of free energy due-to 
diffusion, in the cell, of the second solute is 


CF m0 tit dost COW C*)? ae leD eV Ce ¥C,* 
dt Ci ae dV. (77) 
V V 


C,' 


In (77) the integrand of the first integral is positive, that of the sec- 
ond negative. The second integral is thus the deficit in decrease in 
free energy, per unit time, due to the presence of the first solute 
which is being maintained in a stationary distribution. The second 
integral is thus the “cost of maintenance” (free energy per unit time) 
of the non-equilibrium concentration distribution of the first solute. 
Thus one may say that (73) gives the “cost” of concentrating the 
first solute in the cell and the second integral in (77) gives the “cost” 
per unit time thereafter to maintain the configuration of C,' in the 
cell. It is not necessary that q, = 0 in the above but this simple case 
is clearest. It will be shown (1950b) that a similar situation holds 
for a cell producing or consuming any number of substances at any 
constant rates g,;. Cases will also be presented in which the q, are 
neither constant nor independent. The above discussed situation is 
not offered as a “formal explanation” of anything. It does have these 
striking, if obvious, features: One meets the classical situation of a 
solute to which the membrane is freely permeable existing in higher 
concentration in the cell than in the environment. This is a non-equi- 
librium situation the existence and maintenance of which depends 
upon the metabolic state of the cell (substances being metabolized 
and the rates). If the metabolism of the cell is altered the non-equi- 
librium configuration is altered; in the extreme it vanishes if me- 
tabolism is stopped. This situation is a strict consequence of the laws 
of diffusion and requires no “one-sided” permeabilities, electrostatic 
forces etc. for explanation. It is not implied here that all such situa- 
tions result from the laws of diffusion alone but the relevancy of the 
above remarks is plain. 

Since the gradients in the above case are due to spontaneous 
chemical reaction in the cell one may by specifying the reactions com- 
pute under simple assumptions the efficiency of the use of “metabolic 
energy” for transport and concentration but this cannot be discussed 


here. 
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5. General Remarks. 


It is not proposed to discuss here in detail the question of 
whether thermodynamics in the usual sense of the word is applicable 
to non-equilibrium systems and irreversible processes but a few re- 
marks seem necessary. It is indeed true that thermodynamics is most 
frequently applied to such systems. Every calculation of 4F from 
the familiar equation (12) of III is such a case. The term “equilib- 
rium” as frequently used in thermodynamics does not imply a state 
such that the system left to itself will not change but a state describ- 
able in terms of thermodynamic or macroscopic variables (cf. discus- 
sion of LaMer et al, 1949). This may necessitate thermal and me- 
chanical (but not chemical) equilibrium in the true sense of the word. 
For the system discussed in III there was assumed an assignable 
value of T and P for the whole system; if the mole numbers 7” can 
be prescribed the description is complete. If the mole numbers change 
in such a way that they can be specified at each stage of the process, 
the change in free energy can be calculated for any change and the 
time rate of change of the free energy (or any other function of 
state), cf. (1) and (6) of III, can be computed and there is nothing 
else to discuss. If T, P and C; are not uniform in the system then the 
system must be divisible into parts in which these variables are uni- 
form; these parts may be infinite in number, cf. (47), requiring a 
space integral for the complete description. Each element is in “‘equi- 
librium” in the sense that there are assignable values of 7, P and C,, 
which obtain throughout the element. In extreme cases an adequate 
description may not be possible even by statistical methods but these 
are unusual systems. This latter point is well known and has been 
discussed by P. W. Bridgman (1941). The difficulties when they do 
arise are in general not associated with the fact that the system is 
changing at a finite rate but with such matters as the choice of suit- 
able macroscopic variables and the difficulty of an operational defini- 
tion of temperature as discussed by Tolman and Fine (1948). J. M. 
Reiner and S. Spiegelman (1945) in discussing Eckart’s (1940b) 
equations give a completely misleading discussion (p. 81-82) imply- 
ing that the formulations of thermodynamics are not available for 
irreversible processes. Their point of view is rather widespread. They 
further state that only inequalities are available from thermodynam- 
ics for the discussion of irreversible processes. That this is not cor- 
rect is abundantly clear from section 1 of this paper, equation (6), 
and from Eckart’s own work. They also, as a result of this point of 
view, fail to recognize the conclusions from (21) and (22) of III and 
(43) of this paper, which equivalent conclusions Eckart also proved 
in a different manner, and give an erroneous discussion (p. 91) of 
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entropy increase and energetic coupling. A. Rosenbluth, N. Weiner, 
W. Pitts and J. Ramos (1948) have, most elegantly, stated that ther- 
modynamics can be applied to systems near equilibrium and which 
change slowly but that to do so is an approximation. It appears to 
this writer that the approximation, as discussed above, must be that 
of describing the system and not associated with change per se, rapid 
or otherwise. 

We hope not to have implied that there are no such problems to 
be met but the unfortunate stress of “equilibrium” and reversibility 
in connection with the applications of thermodynamics seems great 
enough to merit discussion. 


The author is deeply indebted to Miss June Kosarko for the prep- 
aration of this manuscript. 
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BOOK REVIEW 


GEORGE BOoLe. The Mathematical Analysis of Logic. 1948. 82 pp. New York: 
Philosophical Library. (First published in 1847. Cambridge: Macmillan, 
Barclay, & Macmillan; London: George Bell.) $3.75. 


The re-publication of George Boole’s classical pamphlet, The Mathematical 
Analysis of Logic 101 years after its first appearance marks a recognition of the 
importance achieved by so-called “symbolic logic” in modern science. 

The various applications of Boolean algebra or “logical calculus” of interest 
to mathematical biology should be familiar to the readers of this journal through 
the 1948 paper of McCulloch and Pitts: and the subsequent development of the 
ideas embodied therein. 

To this reviewer it seems that Boole’s most significant contribution to science 
has been not so much through his invention of a particular logico-symbolic tech- 
nique as through his emphasis on the symbolic nature of both logic and mathe- 
matics. It is this emphasis which brings out the inherent kinship of the two dis- 
ciplines. Boole’s voice was one of the first to declare clearly that it is not fruitful to 
define mathematics as the science of quantity or magnitude; that although mathe- 
matics does deal most effectively with magnitude, it need not be its sole subject 
matter. Today we would define the subject matter of mathematics as structure, 
and such a notion, although not yet explicitly stated, seems to inspire Boole’s ap- 
proach. 

To those who are already familiar with the simple scheme of Boolean algebra, 
the most interesting part of Boole’s book will be the Introduction. Two theses are 
defended there, (1) the fundamental unity of science and (2) the close relation 
between logic and mathematics. The implication of this orientation can be felt, 
namely, that the logicomathematical method is destined to be language of unified 
science. By a strange coincidence some of the most vociferous objections to this 
orientation came from Sir W. Hamilton, a logician bearing the same name as the 
scientist whose contributions to mathematical physics were enormous. Some older 
arguments against mathematicising knowledge are also taken up and disposed of. 
All of them reveal the ancient dichotomy between “how” and “why,” which for- 
tunately is now fairly thoroughly eliminated from science but unfortunately still 
lurks in undisinfected corners of philosophy. 

One cannot, of course, reproach Boole, unaware as he was of modern seman- 
ties, that he failed to dispose of the dichotomy between “how” and “why” by sub- 
jecting it to operational analysis. Instead he did something else. He accepted the 
verdict of those philosophers who declared that since philosophy is concerned with 
“real existence” and the “research of causes,” matters presumably too lofty to 
include logical analysis, therefore logical analysis forms no part of philosophy. 
In view of this verdict, George Boole concluded that it was justified for logic, 
exiled from the Elysian fields of philosophy, to seek refuge in mathematics. 

And a very good refuge it turned out to be. 
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